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B.  ABSTRACT 
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ONR  Grant  ONR  N00014-07-1 -0597  has  supported  a  four  phase  program.  Under  Phase  (1),  there  were  three 
materials  thrusts,  (a)  The  first  was  on  the  continued  development  of  hybrid  pulse  laser  deposited  and  liquid  phase 
epitaxy  films  of  low  loss  cubic  yttrium  iron  garnet  (YIG)  and  lithium  ferrite,  uniaxial  barium  hexaferrite,  and  planar 
zinc-Y  hexaferrite.  (b)  The  second  was  on  a  renewed  effort  on  bulk  and  thick  film  ultra-dense  fine  grain 
poly  crystalline  YIG  and  spinel  ferrite  materials,  and  on  native  oxide  multilayers,  (c)  The  third  was  on  high 
magnetization  metallic  films,  notably  of  Permalloy,  iron,  and  FeTiN.  All  of  these  classes  of  materials  hold  promise 
for  different  defense  related  radar  and  signal  processing  device  applications.  The  ultimate  applications  include  high 
Q  filters,  phase  locked  microwave  pulse  sources,  microwave  and  millimeter  wave  devices  such  as  isolators, 
circulators,  phase  shifters,  secure  signal  processing  devices,  and  wide  band  tunable  filters.  Under  Phase  (2),  the 
focus  was  on  magnetic  loss  properties,  both  on  and  away  from  the  peak  in  the  ferromagnetic  resonance  (FMR) 
response,  and  the  change  in  the  microwave  response  at  high  power.  These  studies  are  important  because  the  low 
power  losses  control  device  insertion  loss,  while  the  high  power  response  sets  the  limit  on  input  power.  These 
properties  control  the  viability  of  a  given  material  for  useflil  devices.  Under  Phase  (3),  the  Magnetics  Laboratory  at 
Colorado  State  University  (CSU)  worked  with  other  colleagues  and  collaborators  in  sister  ONR  programs  at  the 
University  of  California  -  Berkeley,  Georgia  Institute  of  Technology,  and  Oakland  University  to  combine  the 
magnetic  systems  developed  and  characterized  at  CSU  with  ferroelectric  films  and  systems  to  achieve  new 
multifunctional  properties.  Under  Phase  (4),  the  unique  expertise  of  the  CSU  facility  in  Brillouin  light  scattering 
was  applied  to  the  BiFe03  system,  among  others,  to  study  the  basic  magnetic  excitation  properties  of  these  new  and 
exciting  multi-ferroic  systems  that  exhibit  both  ferroelectricity  and  canted  ferrimagnetic,  or  weak  ferromagnetic, 
order.  One  can  envision  the  use  of  electric  fields  to  tune  the  spin  wave  response  in  this  system  and  produce  a  new 
class  of  low  power  consumption  tunable  microwave  and  millimeter  wave  devices. 
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C.  RESULTS 


1.  Overview 

There  are  few  centers  of  expertise  in  the  world 
which  are  capable  of  quality  research  and  development 
work  in  the  area  of  microwave  magnetic  materials.  The 
Magnetics  Laboratory  in  the  Department  of  Physics  at 
Colorado  State  University  (CSU),  Fort  Collins,  is  well 
equipped  for  a  wide  range  of  high  frequency  magnetic 
measurements  at  low  and  high  power,  at  low  and  high 
field,  over  a  wide  range  of  frequencies,  and  as  a 
function  of  temperature.  This  team  also  has  the  resident 
expertise  to  select  and  investigate  the  critical  materials 
problems  that  are  relevant  to  the  needs  cited  above 
Much  of  this  infrastructure  has  been  established  over 
the  past  decade  with  Office  of  Naval  Research  (ONR) 
support. 

Past  CSU  work  has  addressed  and  solved  numerous 
microwave  loss  and  materials  problems  as  they  relate  to 
both  fundamental  understanding  and  device  needs.  In 
the  area  of  ferrites,  these  include  (1)  the  role  of 
microstructure  in  the  low  and  high  power  loss 
properties  of  polycrystalline  yttrium  iron  garnet  (YIG), 
(2)  the  low  and  high  power  microwave  properties  of 
substituted  lithium  ferrite  materials,  (3)  the  microwave 
properties  of  arc  plasma  spray  lithium  ferrite,  (4)  the 
origins  of  the  large  losses  in  hexagonal  ferrite  materials 
for  millimeter  wave  applications,  (5)  the  high  power 
microwave  properties  of  hexagonal  ferrite  materials,  (6) 
the  characterization  of  liquid  phase  epitaxy  yttrium  iron 
garnet  (YIG)  films  produced  with  special  fluxes,  (7)  the 
characterization  of  ultra  dense  polycrystalline  ferrites 
for  microwave  applications  prepared  by  hot  isostatic 
pressing  (HIPPING)  techniques,  (8)  the  growth  and 
characterization  of  pulsed  laser  deposited  (PLD)  ferrite 
films  of  YIG  and  barium  hexaferrite  with  losses  which 
are  as  good  as  the  best  bulk  single  crystals,  and  (9)  the 
successful  PLD  growth  of  low  loss  zinc  lithium  ferrite 
films. 

In  the  area  of  metallic  ferromagnetic  films,  which 
represents  a  promising  system  for  wide  band  tunable 
filters  in  the  microwave  and  millimeter  regime,  for 
example,  CSU  has  been  a  key  contributor  to  (1)  the 
understanding  of  phenomenological  damping  in  metal 
films,  (2)  the  elucidation  of  microwave  loss  properties 
and  most  recently,  (3)  a  new  understanding  of  the  high 
power  properties  as  well.  A  significant  part  of  this 
work  has  been  accomplished  during  the  past  and  current 
grant  periods  under  ONR  support. 


In  addition  to  the  above  microwave/millimeter 
wave  materials  characterization  and  device  physics 
development  work,  the  Colorado  State  University 
program  has  provided  numerous  contributions  that  have 
advanced  the  understanding  of  microwave  loss 
processes  and  nonlinear  spin  wave  interactions,  both  for 
ferrites  in  general  and  for  thin  films  in  particular. 
These  include  (1)  a  general  formulation  of  the  theory  of 
magnetostatic  waves  in  anisotropic  magnetic  materials, 
(2)  theoretical  analysis  of  spin  wave  instability 
processes,  both  first  and  second  order,  for  materials 
with  a  general  ellipsoidal  shape,  a  general  anisotropy, 
and  a  general  pumping  field  configuration,  (3)  direct 
identification  of  the  spin  wave  interactions  responsible 
for  the  onset  of  nonlinear  loss  in  ferrites  at  high  power, 
and  (4)  practical  theoretical  models  of  the  two  magnon 
scattering  interaction  and  calculations  of  the  resulting 
linewidths  and  off  resonance  losses  in  ferrite  materials. 
A  significant  part  of  this  work  has  also  been 
accomplished  under  ONR  support. 

The  previous  phase  of  this  continuing  ONR 
program  under  N000 14-06- 1-0889  reported  new  results 
on  the  static  and  microwave  properties  of  composite 
ferrite-ferroelectric  materials,  nonlinear  and  high  power 
properties  of  ferrite  and  metallic  films,  low  field  loss 
properties  of  microwave  ferrites,  and  new  theoretical 
approaches  for  the  analysis  of  microwave  loss  in 
magnetic  materials.  See  Final  Report  dated  29 
November  2007. 

This  phase  of  the  ONR  program  achieved  new 
results  in  several  areas.  The  results  include  the 
experimental  and  theoretical  verification  of  loss 
processes  in  metallic  films,  the  development  of  in-plane 
c-axis  barium  hexaferrite  films  for  millimeter  wave 
device  applications,  the  demonstration  of  a  record 
electric  field  tuning  response  for  the  ferromagnetic 
resonance  line  in  a  ferrite  yttrium  iron  garnet  (YIG) 
film  -  ferroelectric  barium  strontium  titanate  film 
hetero structure,  new  results  on  loss  mechanisms  in  utra- 
dense  bulk  polycrystalline  YIG  materials,  and  the  first 
experimental  confirmation  of  nonlinear  three  magnon 
confluence  processes  in  ferrite  films.  The  presentation 
of  results  in  this  Final  Report  is  intended  to  be  brief  and 
succinct.  Section  C.2  lists  personnel  supported,  in 
whole  or  in  part,  by  category.  Section  C.3  contains  a 
list  of  the  six  archival  publications  realized  during  this 
funding  period. 
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2.  Education  and  Human  Resources 

Personnel  supported  in  whole  or  in  part  and 
degrees  granted  during  the  current  grant  period  are 
indicated  below: 

Visiting  scientists:  1 

Senior  Staff  2 


Postdoctoral  fellows:  4 

Undergraduate  students:  3 

High  school  summer  apprenticeships:  6 

Names,  degree  specifics,  and  dates  are  available  on 

request. 
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PHYSICAL  REVIEW  B  77,  054427  (2008) 


Microwave  damping  in  polycrystalline  Fe-Ti-N  films:  Physical  mechanisms  and  correlations  with 

composition  and  structure 

Sangita  S  Kalarickal,1  Pavol  Krivosik,1 2  Jaydip  Das,1  Kyoung  Suk  Kim,3  and  Carl  E.  Patton1 
1  Department  of  Physics,  Colorado  State  University,  Fort  Collins,  Colorado  80523 ,  USA 
2Slovak  University  of  Technology,  81219  Bratislava,  Slovakia 
}Deparment  of  Materials  Science  ami  Engineering,  Korea  University,  Seoul  136-701,  Korea 
(Received  30  October  2007;  published  21  February  2008) 

Ferromagnetic  resonance  (FMR)  derivative  linewidths  were  measured  from  3  to  12  GHz  on  50  nm  thick 
sputtered  polycrystallinc  Fe-Ti-N  films  with  3  at.  %  titanium  and  a  nitrogen  content  (jc^)  from 
1.9  to  12.7  at.  %.  The  measurements  were  made  with  both  stripline  and  waveguide  FMR  spectrometers  Line- 
widths  were  generally  lowest  at.*v=7  at.  %,  with  derivative  linewidth  (A H)  values  in  the  15-25  Oe  range  and 
a  nominally  linear  increase  with  frequency  (/).  This  minimum  linewidth  composition  is  connected  with  the  bcc 
to  bet  structural  transition  in  the  Fe-Ti-N  system.  Linewidths  increased  at  both  larger  and  smaller  jc^  values  and 
were  accompanied  by  the  development  of  a  more  rounded  frequency  profile  that  is  indicative  of  two-magnon 
scattering.  All  of  the  A H  vs  /  data  could  be  fitted  successfully  with  a  constant  inhomogcncity  hroadcning 
linewidth  of  8-1 1  Oc,  a  two-magnon  scattering  (TMS)  linewidth  from  the  random  grain-to-grain  fluctuations 
in  the  effective  anisotropy  held  directions  for  the  polycrystal,  and  a  magnon-electron  (m-e)  intrinsic  relaxation 
term  modeled  through  Gilbert  damping  with  a  single  a  value  of  0.003.  The  actual  hts  were  done  through  the 
convolution  of  a  Gaussian  linewidth  for  the  inhomogeneity  term  and  a  Lorentzian  linewidth  for  the  TMS  and 
m-e  terms.  The  fitted  anisotropy  field  parameters  from  the  TMS  analysis  ranged  between  398  and  883  Oe,  with 
the  minimum  also  at  the  bcc  to  bet  structural  transition  at  %=7  at.  %. 

DOI:  10. 1 1 03/PhysRevB. 77.054427  PACS  numbers):  76.50 -fg,  75,50.Bb,  75.30,Ds 


I.  INTRODUCTION 

Fe-X-N  thin  films,  where  X  denotes  Al,  Ti,  or  Ta,  for 
example,  have  been  a  subject  of  recent  interest  due  to  the 
large  saturation  magnetization,  high  initial  permeability,  and 
low  coercive  force  in  this  system.1-8  Nanocrystalline  Fe-X-N 
alloy  films,  moreover,  provide  an  attractive  option  for  high 
frequency  and  high  speed  applications  for  microwave  de¬ 
vices  and  information  storage.  While  some  important  con¬ 
nections  between  film  composition,  crystal  structure,  micro¬ 
structure,  static  magnetic  properties,  and  magnetodynamic 
properties  have  been  reported,2*3,7-9  there  are  serious  gaps  in 
the  understanding  of  the  microwave  loss. 

The  goal  of  the  present  work  was  to  study  the  microwave 
damping  properties  in  soft,  polycrystalline  Fe-Ti-N  thin  films 
as  a  function  of  frequency  and  nitrogen  content.  This  work 
connects  the  change  in  structure  and  anisotropy  with  nitro¬ 
gen  to  the  magnetodynamics  and,  in  particular,  to  the  ferro¬ 
magnetic  resonance  (FMR)  linewidth.  Specifically,  the  room 
temperature  FMR  response  was  measured  as  a  function  of 
frequency  from  3  to  12  GHz  for  in-plane  magnetized  films 
with  3  at.  %  titanium  and  1.9-12.7  at.  %  nitrogen.  The  line 
shapes  show  a  more  Gaussian  than  Lorentzian  character  at 
low  frequency  and  a  more  Lorentzian  than  Gaussian  charac¬ 
ter  at  high  frequency.  This  observation  leads  to  the  use  of  a 
convolution  approach  to  the  overall  linewidth  analysis.  The 
linewidth  data  were  fitted  to  a  combination  of  three  pro¬ 
cesses:  line  broadening  induced  by  inhomogeneities,  two- 
magnon  scattering  due  to  randomly  oriented  grains  in  the 
polycrystal,  and  intrinsic  Gilbert  damping.  The  analysis  indi¬ 
cates  that  these  processes  quantitatively  define  the  overall 
loss  picture  and  that  two-magnon  processes  dominate.  The 
extracted  fit  parameters  correlate  well  with  the  expected 


structural  changes  with  nitrogen  content  and  static  magnetic 
parameters. 

The  paper  is  organized  as  follows.  Section  II  gives  a  brief 
review  of  past  work  on  the  Fe-Ti-N  thin  film  system  Section 
III  describes  the  thin  film  materials  and  gives  details  on  the 
experimental  setup.  Section  IV  presents  the  experimental 
FMR  linewidth  results.  This  includes  example  absorption  de¬ 
rivative  profiles  and  Gaussian  and  Lorentzian  best  fits  to 
these  profiles,  as  well  as  data  on  linewidth  vs  frequency  and 
nitrogen  content.  Section  V  considers  line  shapes  and  line- 
widths  for  inhomogeneous  line  broadening  and  bona  fide  re¬ 
laxation  processes,  reviews  the  theoretical  basis  for  the 
mechanisms  of  choice  used  here,  and  gives  working  equa¬ 
tions  for  the  linewidth  analysis.  Section  VI  presents  theoret¬ 
ical  fits  to  the  data  and  discusses  the  different  linewidth  con¬ 
tributions  as  they  relate  to  inhomogeneities,  two-magnon 
scattering,  and  intrinsic  losses.  Section  VII  provides  a  sum¬ 
mary  and  conclusion. 

II.  BRIEF  NOTE  ON  STRUCTURE  AND  LINEWIDTH 
CONNECTIONS 

Previous  work  on  Fe-X-N  alloy  films  has  shown  that  there 
is  a  transition  from  a  body  centered  cubic  (bcc)  structure  to  a 
body  centered  tetragonal  (bet)  structure  at  a  nitrogen  concen¬ 
tration  (jty)  of  about  7  at.  %.3*4, 10,11  jn  a  recent  paper,  Das  et 
al.12  have  also  established  connections  between  the  structural 
transition  and  a  concomitant  decrease  in  cubic  anisotropy 
and  increase  in  uniaxial  anisotropy  as  the  nitrogen  level  is 
increased.  The  operational  anisotropy,  cubic  or  uniaxial,  is 
minimum  at  around  7  at.  %. 

While  the  above  results,  especially  the  data  from  Ref.  12, 
demonstrate  the  structural,  crystallographic,  and  static  mag- 
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TABLE  1.  Summary  of  parameters  for  the  Fe-Ti-N  films. 


Nitrogen 

content 

xN 

(at.  %) 

Grain 

size 

f 

(nm) 

Saturation 

induction 

4ttA/v 

(kG) 

FMR 
linewidth 
at  9.5  GHz 

A  H 

(Oe) 

19 

20 

19.0 

76 

3.9 

15 

18.8 

57 

5.4 

10 

18.2 

46 

7 

8.5 

16.8 

27 

8.4 

7.5 

15.3 

27 

12.7 

4 

13.9 

32 

netic  property  connections,  similar  links  with  dynamic  prop¬ 
erties  are  not  yet  clear.  These  systems  do  show  relatively  low 
FMR  linewidths9  that  (1)  decrease  as  xN  increases  and  (2) 
appear  to  correlate  with  the  corresponding  decrease  in  grain 
size.  References  8  and  1 2  also  suggest  a  possible  connection 
between  linewidth  and  magnetization  ripple  in  thin  films. 
The  above  cited  FMR  linewidth  work,  although  somewhat 
limited,  provides  a  natural  lead  into  the  present  work.  These 
results,  described  below,  elucidate  the  connections  between 
grain  size,  structural  changes,  and  linewidth.  To  a  large  ex¬ 
tent,  these  results  also  resolve  the  predominant  two-magnon 
and  inhomogeneity  broadening  contributions  to  the  FMR 
linewidth  in  the  Fe-Ti-N  thin  film  system. 

III.  SAMPLE  PROPERTIES  AND  FERROMAGNETIC 
RESONANCE  MEASUREMENT  TECHNIQUES 

As  with  the  work  in  Ref.  12,  the  polycrystalline  Fe-Ti-N 
films  were  provided  by  C.  Alexander,  Jr.  of  the  University  of 
Alabama.  The  50  nm  thick  films  were  prepared  by  dc  mag¬ 
netron  sputtering  on  glass  substrates.  A  dc  magnetic  field  of 
about  300  Oe  was  applied  in  the  plane  of  the  substrates  dur¬ 
ing  sputtering  to  induce  uniaxial  anisotropy.  All  films  had 
about  3  at.  %  Ti,  The  samples  had  xN  values  of  1.9,  3.9,  5.4, 
7.0,  8.4,  and  12.7  at.  %.  The  compositions  were  determined 
by  x-ray  photoelectron  spectroscopy.  The  average  grain  size, 
as  detenu ined  by  transmission  electron  microscopy,  ranged 
from  20  nm  at  xN=  1 .9  at.  %  to  about  4  nm  at  xN 
=  12.7  at.  %.  Preparation  details  are  given  in  Ref.  2.  A  thor¬ 
ough  discussion  of  the  static  magnetic  properties  and  con¬ 
nections  with  structure  is  given  in  Ref.  12. 

Table  I  lists  key  sample  parameters.  The  xN  and  grain  size 
(£)  values  are  from  Ref.  2.  The  saturation  induction  (47tA/v) 
values  are  from  Ref.  12.  The  listed  9.5  GHz  FMR  derivative 
linewidth  (A H)  values  are  from  the  present  measurements 
Note  that  saturation  induction,  grain  size,  and  linewidth  val¬ 
ues  generally  decrease  with  increasing  nitrogen  content.  This 
listing  also  demonstrates  the  relatively  low  A H  values  for  the 
Fe-Ti-N  film  system. 

Two  systems  were  used  for  the  linewidth  measurements,  a 
broadband  stripline  transmission  spectrometer13,14  for 
2-7  GHz  and  a  shorted  waveguide  reflection 
spectrometer15,16  for  8-12  GHz.  The  FMR  absorption  de- 


R2 

3  GHz 

5  25  GHz 

Lorentzian 

0  9512 

0  9968 

Gaussian: 

0.9847 

0.9600 

q  — ■ — i — » — - — ■ — - — i — * — - — * — - — 1 — 

50  75  100 

Applied  field  (Oe) 


FIG.  1.  (Color  online)  (a)  and  (b)  show  example  ferromagnetic 
resonance  absorption  derivative  vs  field  profiles  for  the  film  with 
^=5.4  at.  %  at  3.0  and  5.25  GHz,  respectively.  The  dashed  (red) 
and  solid  (blue)  curves  show  Gaussian  and  Lorentzian  fits  to  the 
data.  The  top  table  lists  the  regression  parameters  (R  )  for  the  four 
fits. 

rivative  vs  field  profiles  were  generally  undistorted.  The 
FMR  derivative  linewidth,  taken  here  as  A 7/,  was  measured 
as  the  difference  between  the  extrema  of  the  FMR  profile. 
The  variation  in  these  linewidths  with  frequency  and  nitro¬ 
gen  content,  along  with  the  analysis  of  these  data  in  terms  of 
basic  loss  and  line  broadening  processes,  is  the  focus  of  this 
work.  It  is  to  be  emphasized  that  all  experimental  and  theo¬ 
retical  linewidths  presented  here  are  given  as  derivative  line- 
widths. 


IV.  FERROMAGNETIC  RESONANCE  LINEWIDTH 
RESULTS 

Figure  1  shows  example  FMR  absorption  derivative  pro¬ 
files  for  the  sample  with  ^=5.4  at.  %  at  3.0  and  5.25  GHz. 
The  open  circles  show  the  data.  The  dashed  and  solid  curves, 
respectively,  show  Gaussian  and  Lorentzian  derivative  fits  to 
the  data.  For  this  presentation,  the  standard  regression  pa¬ 
rameter  R  is  used  as  a  measure  of  the  goodness  of  the  fits.  A 
perfect  fit  would  correspond  to  R2=\.  These  profiles  demon¬ 
strate  two  things.  First,  one  can  see  that  the  resonances  are 
reasonably  sharp  and  symmetric,  with  relatively  narrow  line- 
widths,  as  indicated.  The  second  point  concerns  line  shapes. 
One  can  see  that  for  the  3.0  GHz  data  in  Fig.  1  (a),  the  Gauss¬ 
ian  fit  is  better  than  the  Lorentzian  fit,  and  for  the  5.25  GHz 
data  in  Fig.  1(b),  the  Lorentzian  fit  is  better.  While  these 
differences  are  visible  from  the  fits  shown,  the  variation  in 
the  goodness  of  the  fits  shows  up  quite  strongly  in  the  listed 
R2  values.  These  clear  differences  provide  the  motivation  for 
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FIG.  2.  (Color  online)  FMR  derivative  linewidth  A H  vs  fre¬ 
quency  (w/2tt)  for  all  samples.  The  symbols  show  the  data  for 
different  nitrogen  contents,  as  indicated.  The  curves  show  model  fits 
to  the  data  based  on  the  theoretical  considerations  discussed  in  the 
text. 

the  convolution  approach  to  the  linewidth  analysis  used  be¬ 
low. 

Figure  2  shows  the  FMR  derivative  linewidth  A//  as  a 
function  of  frequency  ((dI2tt)  for  all  samples.  The  symbols 
show  the  data  for  the  different  nitrogen  levels,  as  indicated. 
The  curves  show  fits  to  the  data,  based  on  theoretical  con¬ 
siderations  of  two-magnon  scattering,  inhomogeneous  line 
broadening,  and  intrinsic  magnon-electron  (m-e)  scattering 
processes.  The  theory  and  fits  are  discussed  in  Secs.  V  and 
VI.  Note  that  the  frequency  cu,  specified  in  rad/s,  plays  an 
important  role  in  the  theoretical  discussion  to  follow.  The 
frequency  parameter  used  for  presentations  of  the  data  will 
usually  be  given  in  practical  units  as  GHz. 

There  are  two  important  points  to  be  noted  from  Fig.  2. 
First,  it  is  clear  that  the  linewidth  is  not  a  linear  function  of 
frequency.  While  it  is  true  that  single  crystal  films  and  thin 
polycrystalline  films  can  show  a  linear  A//  vs  to 
signature,14-17"19  this  cannot  be  taken  as  a  general 
statement.20"22  From  the  data,  one  can  see  that  the  response 
is  closest  to  linear  for  an  intermediate  nitrogen  level  at  xN 
~7  at.  %  This  composition  also  corresponds  to  the  point 
where  the  linewidth  is  a  minimum  As  one  goes  below  or 
above  this  level,  one  obtains  larger  linewidths  and  a  more 
accentuated  curved  response.  These  effects  are  larger  for  xN 
values  below  7  at.  %  than  for  xN>l  at.  %.  Figure  3  shows 
the  data  of  Fig.  2  in  a  A H  vs  xN  format  for  three  selected 
frequencies,  as  indicated.  In  this  format,  the  data  demonstrate 


FIG.  3.  (Color  online)  FMR  derivative  linewidth  A H  vs  nitrogen 
content  xN  at  three  selected  frequencies,  as  indicated. 
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explicitly  the  linewidth  minimum  at  xN**l  at.  %  and  the 
more  rapid  rate  of  change  for  *iV<7  at.  %  noted  above. 

The  data  shown  in  Figs.  2  and  3  can  be  understood  in 
terms  of  the  static  magnetic  and  structural  properties  of  the 
films.  Inhomogeneous  line  broadening  contributes  a  fre¬ 
quency  and  composition  independent  background  linewidth. 
The  linewidth  minimum  at  xN~l  at.  %  is  related  to  the  es¬ 
tablished  structural  transition  at  this  composition*1  and  a 
minimum  value  for  the  overall  magnetic  anisotropy  of  the 
films.12  The  increases  in  linewidth  above  and  below  the  xN 
~1  at.  %  point  as  well  as  the  enhanced  curvature  of  the  A// 
vs  a)  response  are  due  to  two-magnon  scattering  (TMS).23  26 
These  changes  are  related  to  the  increase  in  anisotropy12  and 
grain  size  changes.  The  remaining  linewidth  component  fol¬ 
lows  the  well  established  Gilbert  (G)  phenomenological 
model  with  a  single,  composition  independent  a  value  of 
0.003.  This  a  value  is  consistent  with  intrinsic  magnon- 
electron  scattering,  the  main  intrinsic  relaxation  process  for 
ferromagnetic  resonance  in  metallic  films.27  30 

The  next  section  provides  an  overview  of  these  three  line- 
width  mechanisms  and  develops  working  equations  for  fits  to 
the  data.  Section  VI  presents  the  results  of  these  fits  and 
discusses  the  role  of  the  different  mechanisms. 

V.  LINEWIDTH  MECHANISMS  AND  WORKING 
EQUATIONS 

The  most  important  point  to  note  at  the  outset  of  this 
section  is  that  inhomogeneities  and  physical  loss  processes 
lead  to  different  line  shapes,  and  one  cannot  simply  add  these 
two  different  types  of  linewidths  to  obtain  a  total  linewidth. 
To  lowest  order,  the  effect  of  inhomogeneities  amounts  to  the 
superposition  of  many  local  FMR  absorption  profiles  with 
some  spread  in  FMR  fields  that  define  the  extent  of  the  in- 
homogeneity.  For  random  inhomogeneities,  the  spread  will 
follow  a  Gaussian  distribution.  The  corresponding  FMR  line- 
width  will,  in  the  small  signal  limit,  also  have  a  Gaussian 
shape. 

The  other  two  processes,  two-magnon  scattering  and  in¬ 
trinsic  magnon-electron  scattering,  comprise  physical  relax¬ 
ation  mechanisms.  In  the  small  signal  limit  and  to  lowest 
order,  these  processes  will  have  Lorentzian  FMR  absorption 
profiles.  Because  of  this  basic  difference,  any  measured  line- 
width  will  be  a  convolution  of  the  Gaussian  profile  associ¬ 
ated  with  the  inhomogeneous  line  broadening  and  the 
Lorentzian  profile  associated  with  the  relaxation  processes. 
This  type  of  convolution  applies,  of  course,  only  if  these 
effects  are  independent  of  each  other.  The  fits  shown  in  Fig. 
1  support  the  convolution  argument.  From  the  linewidth 
analysis  in  Sec.  VI,  it  will  be  shown  that  the  inhomogeneity 
linewidth  contribution  is  larger  than  the  linewidth  associated 
with  actual  relaxation  processes  at  lower  frequencies,  while 
the  opposite  is  true  at  higher  frequencies.  This  is  precisely 
the  indication  from  Fig.  1. 

Based  on  the  above,  consider  two  contributions  to  the 
FMR  linewidth.  There  is  one  linewidth  associated  with  a 
Gaussian  profile  and  related  to  inhomogeneity  (INH)  line 
broadening,  termed  A //1NH.  In  addition,  there  is  a  second 
linewidth,  one  associated  with  a  Lorentzian  (LOR)  profile 
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and  related  to  relaxation  processes,  termed  A HL0R,  In  the 
present  situation,  one  can  write,  by  way  of  approximation, 
A//LOr=A//tk1s+A//c,  where  A/f^s  and  A HG  denote  the 
respective  two-magnon  scattenng  and  magnon-electron  or 
Gilbert  contributions.  Techniques  for  the  convolution  of 
Gaussian  and  Lorentzian  line  shapes  and  the  use  of  Voigt 
functions  are  well  established.31  Here,  however,  the  focus  is 
on  the  linewidths,  not  on  the  overall  shapes  of  the  profiles. 
Working  equations  for  the  approximation  of  the  derivative 
linewidths  may  be  found  in  Ref.  3 1 .  The  combined  deriva¬ 
tive  linewidth  A H  takes  the  form 

+  0-909 A//|NHA//lor  +  0,462A//LOR 
A/^inh  +  0.462 A//lor 

Note  that  Eq.  (1)  is  specific  to  derivative  linewidths.  As  the 
development  below  will  show,  the  convolution  of  the  inho¬ 
mogeneous  and  Lorentzian  contributions  to  the  FMR  re¬ 
sponse  gives  a  good  fit  to  the  overall  data.  Figure  1  has 
already  demonstrated  the  shift  from  Gaussian  to  Lorentzian 
line  shapes  as  a  function  of  frequency.  Linewidth  analyses 
based  on  a  simple  addition  of  A//Lqr  and  A //1NH  are  widely 
used.18,21,32,33  Based  on  the  remarks  above,  however,  it  is 
clear  that  convolution  represents  a  more  consistent  approach 
when  inhomogeneity  and  relaxation  contributions  to  the  line- 
width  are  both  present. 

Now  consider  the  individual  contributions  to  the  overall 
linewidth.  As  noted,  line  broadening  due  to  inhomogeneities 
arises  from  the  simple  superposition  of  several  local  FMR 
profiles  for  different  regions  of  the  sample  that  are  shifted  in 
field  because  of  some  change  in  the  magnetic  properties.  As 
one  example,  noninteracting  regions  with  different  Ms  or  an¬ 
isotropy  values  can  lead  to  such  shifts.  It  is  the  assumed 
random  distribution  in  these  parameters  relative  to  some 
mean  value  that  gives  rise  to  the  overall  Gaussian  character 
of  the  line  shape. 

Inhomogeneous  line  broadening  can  take  many  forms.  In 
polycrystalline  ferrites  with  a  random  distribution  of  single 
crystal  grains,  for  example,  the  line  shape  will  map  the  ran¬ 
dom  orientations  of  the  crystallographic  easy  axes  of  the 
grains.23  To  lowest  order,  the  corresponding  spread  in  FMR 
fields  will  give  a  constant  A//INH  that  is  independent  of  fre¬ 
quency.  This  model  appears  to  apply  to  the  present  data. 
More  complicated  inhomogeneity  line  broadening  models 
have  also  been  proposed.  These  include  a  field  dependent 
line  broadening  effect  related  to  magnetization  ripple.9  A 
brief  discussion  of  the  ripple  option  will  be  given  at  the  end 
of  Sec.  VI. 

In  addition  to  the  conceptually  simple  Gaussian  line 
broadening  effects  discussed  above,  inhomogeneities  can 
also  give  rise  to  a  Lorentzian  linewidth  contribution  related 
to  the  scattering  of  spin  waves  that  are  at  the  same  frequency 
as  the  driven  FMR  mode.  The  relaxation  mechanism  associ¬ 
ated  with  this  process  is  called  two-magnon  scattering.  The 
process  was  discovered  in  the  late  1950s  and  early  1960s  in 
connection  with  FMR  loss  studies  in  bulk  ferrites.13’23,34,35  In 
this  same  time  frame,  TMS  processes  were  also  proposed  as 
a  possible  linewidth  mechanism  in  thin  films.13,36  In  recent 
years,  there  has  been  a  renewed  interest  in  TMS  processes  in 


metallic  ferromagnetic  films,  driven  mainly  by  applications 
in  high  density  and  high  speed  magnetic  recording.2 124, 37-40 

Briefly  stated,  TMS  processes  are  related  to  short  wave¬ 
length  fluctuations  in  the  local  effective  field  in  a  given 
sample  that  break  the  orthogonality  of  the  FMR  mode  and 
other  degenerate  spin  wave  modes.  Examples  can  include  the 
grain-to-grain  fluctuations  in  the  effective  magnetocrystalline 
anisotropy  field  for  a  polycrystalline  sample  with  randomly 
oriented  grains,  small  pores  or  second  phase  material,  and 
surface  defects.  Such  fluctuations  can  lead  to  a  coupling  be¬ 
tween  the  driven  FMR  mode  and  degenerate  spin  wave 
modes,  and  provide  a  channel  for  TMS  relaxation.  The  cor¬ 
responding  linewidth  will  be  strongly  dependent  on  the  den¬ 
sity  of  states  for  degenerate  spin  wave  modes  and,  hence,  the 
shape  of  the  spin  wave  band.  The  strongest  coupling  is  usu¬ 
ally  for  spin  waves  with  wavelengths  that  are  well  above  the 
size  of  the  inhomogeneity.  The  process  cuts  off  for  wave¬ 
lengths  much  shorter  than  the  inhomogeneity  size.  For  the 
5-20  nm  or  so  grain  sizes  applicable  here,  the  cutoff  will  be 
for  a  spin  wave  wave  number  k  on  the  order  of 
106-107  rad/cm.  This  is  precisely  the  range  of  the  available 
modes  for  the  present  Fe-Ti-N  films  at  microwave  frequen¬ 
cies. 

The  thin  film  TMS  theory  used  here  is  based  on  the  clas¬ 
sical  analysis  developed  in  Ref.  24.  References.  24  and  25 
provide  a  qualitative  introduction  to  thin  film  TMS  processes 
in  general,  details  on  the  theory,  and  working  equations  for 
practical  calculations  and  fits  to  linewidth  data.  Brief  com¬ 
ments*  are  given  below.  Working  equations  for  the  linewidth 
calculations  are  given  in  the  Appendix. 

For  a  tractable  analysis,  one  usually  invokes  the  thin  film 
approximation  introduced  by  Harte,41  in  which  the  spin  wave 
modes  are  taken  to  have  no  spatial  variation  over  the  film 
cross  section.  In  this  limit,  the  thin  film  TMS  linewidth  de¬ 
veloped  in  Ref.  25  may  be  cast  in  the  form 

j  \y\H2t2  f 

AW™s  =  irLr7T  (2) 

v3  J 

Here,  HA  denotes  an  effective  anisotropy  field  for  the  ran¬ 
domly  oriented  grains  defined  through  HA  =  2KI where  K 
is  an  as  yet  unspecified  anisotropy  energy  density  for  the 
material  of  interest.  In  Sec.  VI,  both  cubic  and  uniaxial  an¬ 
isotropy  contributions  will  be  considered.  The  £  parameter 
denotes  the  mean  grain  size.  The  PA(o>)  function  is  equal  to 
/d\y\H)  evaluated  at  H=HhMR,  where  ro^R  is*  the 
field  dependent  FMR  angular  frequency  and  //^mr  is*  the 
resonance  field  at  the  operating  frequency  The  y  parameter 
denotes  the  gyromagnetic  ratio.  The  PA  divisor  in  Eq.  (2) 
accounts  for  the  conversion  between  frequency  and  field 
swept  linewidths.14  The  Ck(£)  is  the  Fourier  transform  of  the 
grain-to-grain  dynamic  internal  field  correlation  function  for 
the  sample.  Detailed  considerations  of  Ck(£)  based  on  the 
anisotropy  in  a  system  of  randomly  oriented  grams  are  given 
in  Ref.  25.  This  Ck(£)  provides  a  specific  form  for  the  k 
dependence  of  the  scattering  for  a  given  £  value.  The  Aok 
takes  into  account  the  averaging  of  the  anisotropy  axis  fluc¬ 
tuations  over  the  sample  and  the  ellipticity  of  the  precession 
response  for  both  the  uniform  mode  and  the  spin  wave 
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modes.  The  integral  folds  in  the  scattering  from  the  uniform 
mode  to  all  available  degenerate  spin  wave  modes.  The  delta 
function  serves  to  select  out  these  degenerate 

modes. 

A  full  discussion  of  the  dependences  of  the  two-magnon 
hnewidth  on  frequency,  grain  size,  and  other  magnetic  prop¬ 
erties  is  beyond  the  scope  of  this  paper.  The  key  points  for 
the  present  work  are  the  following:  (1)  the  different  terms  in 
Eq.  (2)  can  lead  to  a  complicated  dependence  of  the  two- 
magnon  linewidth  with  frequency  and  grain  size,  and  (2)  the 
A//Tms(w)  response  is,  in  general,  not  linear.  As  shown  in  the 
next  section,  the  explicit  nonlinear  A//™s(a>)  response  for 
these  Fe-Ti-N  films  plays  a  crucial  role  in  the  explanation  of 
the  overall  linewidth 

While  two-magnon  scattering  is  a  physical  relaxation  pro¬ 
cess,  it  is  based  on  the  presence  of  inhomogeneities  and  can 
be  taken  as  nonintrinsic.  The  last  process  needed  for  the 
present  analysis  is  m-e  scattering.  This  is  the  sole  intrinsic 
relaxation  process  considered  here.  From  the  m-e  consistent 
Gilbert  model,  one  obtains  a  linewidth,  taken  here  as  A HG, 
that  takes  the  form 


A  Hc  = 


2  aa) 

\3  H 


(3) 


The  linear  increase  in  linewidth  with  frequency  implicit  in 
the  Gilbert  form  has  been  widely  observed  in  ferromagnetic 
metals.14-19  As  noted  above,  the  Gilbert  phenomenology  is 
also  consistent  with  the  physical  relaxation  process  of 
magnon-electron  scattering.  Recent  ab  initio  band  structure 
calculations  of  the  m-e  relaxation  rates  for  iron,  cobalt,  and 
nickel  give  Gilbert  damping  a  values  that  are  consistent  with 
the  present  results  as  well  as  with  the  other  data  cited 
above.28  The  \3  divisor  in  Eqs.  (2)  and  (3)  accounts  for  the 
conversion  between  a  half-power  linewidth  and  the  deriva¬ 
tive  linewidth  considered  here,  under  the  assumption  of  a 
Lorentzian  absorption  profile. 


VI.  LINEWIDTH  ANALYSIS  AND  PHYSICAL 
CONNECTIONS 

Section  V  reviewed  the  INH,  TMS,  and  G  contributions  to 
the  FMR  linewidth,  with  an  emphasis  on  the  possible  fre¬ 
quency  dependences  of  these  linewidth  terms  as  they  relate 
to  the  present  data.  This  section  presents  explicit  theoretical 
fits  to  the  linewidth  data  from  Sec.  IV.  The  results  provide  a 
self-consistent  explanation  of  the  FMR  linewidth  for  this  se¬ 
ries  of  Fe-Ti-N  thin  films.  The  overall  fits  are  based  on  the 
convolution  of  the  Gaussian  inhomogeneity  line  broadening 
term  with  a  combined  IvOrentzian  term  comprised  of  two- 
magnon  scattering  and  magnon-electron  relaxation.  The  cor¬ 
relation  between  the  extracted  TMS  anisotropy  scattering 
field  parameter  and  composition  provides  particularly  strong 
evidence  for  the  importance  of  two  magnon  scattering  in 
these  films. 

Figure  4  shows  the  A H  vs  frequency  data  from  Fig.  2  for 
xN=  1 .9  at.  %,  along  with  computed  curves  for  the  INH, 
TMS,  and  G  hnewidth  contributions  as  well  as  the  total  fitted 
linewidth.  The  open  circles  denote  the  data;  the  dotted, 


FIG.  4  (Color  online)  FMR  derivative  linewidth  AH  vs  fre¬ 
quency  {(o/2tt)  for  xN=  1.9  at.  %.  The  open  circles  show  the  data. 
The  dotted  (blue)  dashed  (red),  and  dot-dashed  (green)  lines  show 
the  computed  inhomogcncity,  Gilbert,  and  two-magnon  scaitcring 
linewidth  components,  respectively.  The  solid  (black)  curve  shows 
the  combined  linewidth  according  to  Eq.  (1). 

dashed,  and  dot-dashed  curves  show  the  INH,  TMS,  and  G 
linewidths,  respectively;  and  the  solid  curve  shows  the  total 
fitted  linewidth  It  was  observed  that  fits  with  Gilbert  a  as  a 
free  parameter  yielded  the  value  of  a  for  all  samples  as 
0.003  ±  0.0002.  The  computed  curves  were  obtained  with  the 
applicable  parameter  values  from  Table  1  and,  as  noted 
above,  a  free  electron  | y| / 2tt  value  of  2.8  GHz/kOe  and  a 
Gilbert  a  value  of  0.003.  Following  Ref.  12,  the  anisotropy 
type  for  .r^/=1.9  at.  %,  needed  for  the  TMS  Aok  function, 
was  taken  as  cubic  and  the  exchange  energy  parameter  A  was 
taken  as  3  X  10~6  erg/ cm.  The  A Him  and  HA  were  taken  as 
fitting  parameters.  The  curves  shown  are  for  A//INH=10  Oe 
and  HA  =  779.7  Oe.  The  regression  parameter  R 2  for  the  fit 
shown  is  0.996. 

Apart  from  the  good  theoretical  fits.  Fig.  4  shows  three 
things.  First,  one  can  see  that  the  convolution  combines  the 
three  component  linewidths  to  give  an  overall  AH  vs  fre¬ 
quency  response  with  a  small  upward  curvature  at  low  fre¬ 
quency.  While  this  is  clearly  a  small  effect,  it  could  tum  out 
to  be  significant  if  the  measurements  could  be  extended  to 
lower  frequencies.  This  and  related  points  will  be  considered 
at  the  end  of  the  section. 

The  second  point  is  that  the  dominant  overall  linewidth 
process,  especially  at  the  higher  frequencies,  is  two-magnon 
scattering.  The  data  and  the  fits  show  that  this  linewidth  con¬ 
tribution  has  a  somewhat  complicated  frequency  depen¬ 
dence,  with  a  linear  increase  at  low  frequency,  a  broad  maxi¬ 
mum  of  8-10  GHz,  and  a  small  drop  at  higher  frequencies. 
The  signature  of  two-magnon  scattering  is  very  different 
from  the  simple  linear  response  implicit  in  the  Gilbert  model. 
Rased  on  the  good  fit  shown  here,  it  is  clear  that  there  can  be 
no  justification  for  the  wide  spread  use  of  the  Gilbert  model 
to  “explain'"  FMR  linewidths  in  metallic  films,  often  without 
any  realistic  consideration  of  the  physical  origins  of  these 
linewidths.  The  third  point  is  that  the  fitted  a  value  of  0.003 
is  consistent  with  expected  values  for  m-e  scattering  in  me¬ 
tallic  films.27-29 

Figure  5  shows  all  of  the  AH  vs  frequency  data  from  Fig. 
2,  broken  out  into  separate  graphs  for  all  the  samples,  along 
with  curves  for  the  overall  linewidth  fits.  The  open  circles 
show  the  data  and  the  solid  curves  show  the  fits.  Apart  from 
the  previously  noted  common  value  of  a  (0.003)  used  for  all 
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F1G.  5.  FMR  derivative  linewidth  A H  vs  frequency  lor  all  the 
samples.  The  solid  curves  give  the  calculations  of  the  linewidth 
with  the  inhomogeneity  line  broadening,  the  two-magnon  scatter¬ 
ing,  and  the  Gilbert  damping  contributions  to  the  total  linewidth. 

fits,  the  curves  in  Fig.  5  are  based  on  the  parameters  from 
Table  1  and  the  additional  linewidth  fitting  parameters  given 
in  Table  TI.  As  with  the  data  in  Fig.  4,  a  constant  background 
inhomogeneity  linewidth  A7/INH  has  been  used  for  each 
graph,  according  to  the  values  in  Table  TI.  Note  that  the  ex¬ 
change  stiffness  constant  A  was  taken  at  3  X  10“6  erg /cm  for 
the  films  with  xN^l  at.  %,12  and  reduced  to  2  X  10-6  and  1 
X  10“6  erg/ cm  for  the  samples  with  .\^=8.4  and  12.7  at.%, 
respectively.  In  view  of  the  bcc  to  bet  structure  change  at 
at.%,  this  change  is  reasonable. 

One  can  also  see  from  Table  TI  that  the  fitted  HA  values 
show  a  minimum  at  xN^l  at.%  and  then  increase  as  one 
goes  above  or  below  this  composition.  Keep  in  mind  that  the 


FIG.  6  (Color  online)  Anisotropy  field  parameter  values  as  a 
function  of  the  nitrogen  content  a  v.  The  solid  (hlaek)  circles  and  left 
side  axis  show  the  HA  values  from  Table  IT.  The  error  hars  are  on 
the  order  of  the  size  of  the  points.  The  open  (blue)  circles  and  right 
side  axis  show  the  averaged  anisotropy  (HA)  values  taken  from  Ref. 
12.  The  solid  curve  serves  as  a  guide  for  the  eye  for  the  trend  of  the 
data. 

two-magnon  fits  are  based  on  cubic  anisotropy  working 
equations  for  values  below  7  at.%  and  on  uniaxial  aniso¬ 
tropy  for  values  of  7  at.%  and  above.  This  is  indicated  in  the 
fifth  column  of  Table  II.  These  anisotropy  details  are  signifi¬ 
cant  for  the  discussion  below. 

The  trend  of  the  data  w  ith  film  composition  was  discussed 
in  connection  with  Fig.  3.  Here,  the  focus  is  on  the  trend  of 
the  fits  with  composition  and  the  two-magnon  scattering  im¬ 
plications.  Note  that  all  of  the  fits  are  quite  good  The  regres¬ 
sion  R2  values  listed  in  Table  II  are  all  reasonably  close  to 
unity.  The  two-magnon  part  of  the  linewidth,  mainly  con¬ 
nected  with  the  A H  vs  frequency  dependences  that  deviate 
from  a  simple  linear  behavior,  starts  out  large  at  xN 
=  1.9  at.%,  initially  drops  off  as  xN  increases,  is  minimum  for 
the  xN=l  and  8.4  at.%  samples,  and  then  becomes  slightly 
larger  as  xN  moves  toward  12.7  at.%. 

The  most  important  result  from  Fig.  5  and  Table  11  relates 
to  the  minimum  two-magnon  linewidth  contnbution  for  xN 
=  7  at.%  and  the  correlations  with  anisotropy  and  structural 
change.  Figure  6  captures  these  results  succinctly.  The  solid 
points  show  the  HA  values  as  a  function  of  x(V  from  Table  II. 
The  open  circles  show  companion  effective  anisotropy  data 
from  Ref.  12.  These  latter  data  derive  from  coercive  force 
measurements  on  these  same  samples  and  extractions  of  an 
effective  cubic  or  uniaxial  anisotropy  energy  density  param- 


TABLE  II.  Summary  of  fit  parameters. 


Nitrogen 

content 

XN 

(at.  %) 

Inhomogcncity 

broadening 

AtflNH 

(Oe) 

Exchange 

constant 

A 

(X  10-6  erg/em) 

Anisotropy  field 

Ha 

(Oc) 

Anisotropy 

type 

Fit  regression 
R 2 

1.9 

10.0±  0.8 

3.0 

779.7  ±4.0 

Cubic 

0.996 

3.9 

8.3  ±0.5 

3.0 

745.1  ±4.1 

Cubic 

0.997 

5.4 

1 1 .0±  0.4 

3.0 

882. 6±  8.6 

Cubic 

0.993 

7 

11. 3  ±0.5 

3.0 

398.4  ±24.9 

Uniaxial 

0.939 

8.4 

10.5  ±0.5 

2.0 

489.4  ±  13.7 

Uniaxial 

0.976 

127 

10  4±  0.7 

1.0 

871.6±  16.1 

Uniaxial 

0.982 
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eter,  denoted  here  as  a  generic  ( K ),  The  ( HA )  values  shown  in 
the  figure  derive  from  an  anisotropy  connection  of  the  form 
(HA)-2(K)f  Ms.  Following  Ref.  12,  and  in  the  spirit  of  the 
two-magnon  analysis  here,  the  values  for  xN^l  at.%  corre¬ 
spond  to  cubic  anisotropy  and  those  for  xN>l  at%  to 
uniaxial  anisotropy.  Note  that  the  right  side  ( HA )  scale  in  Fig. 
6  differs  from  the  left  side  HA  scale  by  a  factor  of  30-80. 
This  difference  in  scale  is  discussed  in  Ref.  12,  The  curve 
through  the  ensemble  of  data  points  is  intended  as  a  guide  for 
the  eye. 

Figure  6  makes  explicit  the  minimum  in  the  HA  values 
from  the  two  magnon  fits  at jtv~7  at.%.  It  also  demonstrates 
the  consistency  between  these  data  and  the  averaged  ( HA ) 
results  from  Ref,  12.  These  results,  in  combination  with  the 
good  match  between  experiment  and  theory  in  Fig.  5,  pro¬ 
vide  convincing  arguments  for  (I)  a  strong  composition  de¬ 
pendent  two-magnon  contribution  to  the  linewidth,  (2)  clear 
correlations  with  a  structure  change  at  xN~l  at.%,  and  (3) 
quantitative  fits  that  support  the  use  of  a  bcc  model  below 
and  a  bet  model  above  the  transition  point.  As  one  such 
quantitative  point  of  comparison,  the  cubic  anisotropy  energy 
density  parameter  from  the  TMS  fit  for  1.9  at%  is  ob¬ 
tained  as  590X  103  erg/cm3,  The  literature  value  for  pure 
iron  is  about  480 X  103  erg/cm3.42 

Up  to  now,  there  has  been  no  comment  on  the  role  of 
grain  size  on  the  calculated  TMS  linewidth.  As  evident  from 
Eq.  (2),  the  grain  size  £  enters  the  TMS  analysis  both  as  a 
prefactor  and  as  a  parameter  inside  the  correlation  integral. 
The  grain  size,  however,  was  not  taken  as  a  fitting  parameter. 
All  of  the  fits  were  based  on  the  measured  grain  sizes  from 
Ref,  12,  as  listed  in  Table  I  It  turns  out  that  any  small  change 
in  grain  size  from  these  measured  values  causes  a  significant 
degradation  in  the  best  fit  matchup  with  the  data  This  means, 
in  turn,  that  the  “best  fit”  £  value  matches  up  nearly  perfectly 
with  the  measured  grain  size.  This  lends  additional  credence 
to  the  two  magnon  linewidth  results  presented  above. 

As  noted  in  Sec.  V,  a  field  dependent  inhomogeneity  line 
broadening  mechanism  related  to  magnetization  ripple  has 
also  been  proposed.  Reasonable  fits  to  the  present  data  could 
also  be  obtained  from  this  model,  the  working  equations  in 
Ref.  9,  and  the  averaged  effective  uniaxial  anisotropy  field 
determinations  from  Ref.  12.  With  this  approach,  the  R2  val¬ 
ues  for  fits  based  on  a  ripple  inhomogeneity  line  broadening 
are  about  the  same  as  those  in  Table  II  for  a  constant  A Him. 
From  these  results,  one  could  argue  that  magnetization  ripple 
is  a  reasonable  source  of  the  A//1NH.  Such  a  mechanism  can 
be  justified  here  because  ripple  line  broadening  has  a  very 
weak  dependence  on  the  applied  field  in  the  FMR  regime  of 
fields  and  frequencies  studied  here.  On  the  other  hand,  the 
ripple  model  also  predicts  a  very  rapid  increase  in  the  inho¬ 
mogeneity  linewidth  at  frequencies  below  2  GHz.  Up  to 
now,  there  is  no  evidence  for  such  a  response  at  low  fre¬ 
quency. 


TABLE  111.  Averaging  coefficients  for  the  variation  in  the  an¬ 
isotropy  direction. 


Averaging  coefficients 

Type  of  anisotropy 

C\ 

c2 

O 

c4 

In-plane  uniaxial 

First  order  cubic 

1/8 

29/420 

3/32 

29/420 

0 

3/140 

1/8 

1/105 

with  nitrogen  concentrations  (jcn)  ranging  from 
1.9-  12.7  at.%.  The  data  show  clear  departures  from  the  lin¬ 
ear  linewidth  vs  frequency  response  that  is  usually  taken  as 
the  norm  for  metallic  films.  The  linewidths  generally  show  a 
minimum  at  the  bcc  to  bet  structural  transition  point  at  xN 
at.%. 

A  model  with  three  terms,  a  constant  inhomogeneous 
broadening  linewidth,  a  two-magnon  scattering  linewidth 
term,  and  an  intrinsic  magnon-electron  term  modeled 
through  a  Gilbert  damping  form,  gives  good  fits  to  the  data. 
The  particular  fits  obtained  here  were  based  on  the  convolu¬ 
tion  of  a  Gaussian  form  for  the  inhomogeneous  term  and  a 
Lorentzian  for  the  two  remaining  physical  processes.  The 
fitted  two-magnon  term  gives  values  of  the  effective  aniso¬ 
tropy  field  in  the  grains  that  are  consistent  with  previous 
reports  and  correlate  with  the  structural  transition  noted 
above. 
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APPENDIX 

This  appendix  gives  working  equations  for  the  evaluation 
of  the  terms  in  Eq.  (2)  for  the  two-magnon  scattering  field 
swept  linewidth  A//TMS.  These  equations  are  specific  to  an 
in-plane  magnetized  thin  film  with  an  average  grain  size  £, 
taken  here  to  be  the  same  as  the  correlation  length  for  the 
grain-to-grain  fluctuation  in  the  effective  anisotropy  field.  All 
of  the  applicable  parameters  have  been  defined  in  the  discus¬ 
sion  of  Eq.  (2).  For  readability,  some  of  these  definitions  are 
repeated  below. 

The  Fourier  transform  of  the  grain-to-grain  effective  field 
correlation  function,  Ck(£),  is  taken  to  have  the  form 


VII,  SUMMARY  AND  CONCLUSION 

In  summary,  the  above  sections  have  described 
3-12  GHz  ferromagnetic  resonance  linewidth  measurements 
and  analyses  for  a  series  of  polycrystalline  Fe-Ti-N  films 


ck(£)  = 


1 

[1  +(*£)2F2’ 


(Al) 


where  A:=  |k|  is  the  spin  wave  wave  number.  In  real  space, 
this  corresponds  to  a  point-to-point  effective  field  correlation 
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that  falls  off  as  e^r~r  K  where  r  and  r'  represent  two  sepa¬ 
rate  positions  in  the  film.24 

The  ellipticity  function  A0u  is  given  by 

c2  I  eo  ek\ 

A()k“cl*,0^k+  +C31  +  )+C4.  (A2) 

¥k  \ek  <V 

Here,  e0  and  ek  denote  the  ellipticity  of  the  uniform  mode 
and  a  given  spin  wave  mode  at  wave  vector  k,  respectively. 
These  ellipticity  factors  are  given  by 


<*0  = 


H  +  47rAf 


(A3) 


and 


*k  = 


H  +  Dkl  +  4irMsNk 


H  +  Dk2  +  4 ttMs{  1  -  Nk)$\n2  0k  ’ 


(A4) 


where  D  is  an  exchange  field  parameter  related  to  the  ex¬ 
change  energy  parameter  A  according  to  D=2A/MS.  The  thin 
film  dipole  field  effective  demagnetizing  factor  Nk  is  given 
by  A*=(l -e~kd)lkdy  where  d  is  the  film  thickness41  and  0k 
is  the  angle  between  the  111-plane  static  held  direction  and  the 
in-plane  wave  vector  k  direction.  The  c j,  c2,  c3,  and  c4  pa¬ 
rameters  are  connected  to  averages  of  the  scattering  over  all 
of  the  randomly  oriented  grains.24  The  numerical  values  of 
these  parameters  for  cubic  and  uniaxial  anisotropies  are 
listed  in  Table  III.  Finally,  the  PA(o))  function  in  Eq  (2)  is 
given  by 


PaM 


(A5) 
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In-plane  c-axis  oriented  barium  ferrite  films  with  self-bias  and  low 
microwave  loss 
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Hybrid  pulse  laser  deposition  and  liquid  phase  epitaxy  methods  have  been  used  to  produce  in-plane 
c-axis  (IPCA)  oriented  barium  ferrite  (BaM)  films  on  a-plane  (1120)  sapphire  substrates  with  low 
microwave  loss  and  a  high  remanence.  Total  thicknesses  were  from  5  to  20  yum  A  reasonable 
compromise  for  low  loss  and  high  remanence  was  reached  at  a  thickness  of  7  yum.  with  a  remanence 
ratio  of  0.84  and  a  59  GHz  peak-to-peak  derivative  linewidth  of  250  Oe.  The  20  /tm  thick  film  had 
a  linewidth  of  110  Oe,  one  of  the  smallest  values  ever  obtained  for  IPCA  BaM  films.  ©  2008 
American  Institute  of  Physics.  [D01‘  10.1063/1.3010374] 


One  important  challenge  for  millimeter  wave  magnetic 
film  device  applications  is  in  the  development  of  low  loss 
self-biased  high  anisotropy  hexagonal  ferrite  films.  Such 
films  have  potential  uses  for  a  wide  range  ot  monolithic  mi¬ 
crowave  integrated  circuit  (MM1C)  magnetic  devices.1-4  The 
critical  parameter  for  self-bias  is  a  high  magnetic  remanence, 
that  is,  a  large  average  magnetization  for  the  film  at  zero 
applied  field.  Low  microwave  loss  is  also  critical  for  a  real 
device  with  an  acceptably  low  insertion  loss. 

Hexagonal  M  type  barium  ferrite  (BaM)  films,  with  a 
high  uniaxial  anisotropy  and  an  effective  anisotropy  field 
(Ha)  of  17  kOe  or  so  and  a  nominal  saturation  induction 
(4vMs)  of  4  kG,  are  good  candidates  for  such  devices.  The 
hexagonal  c-axis  serves  as  the  easy  magnetization  direction. 
The  high  anisotropy  field  allows  one  to  achieve  the  desirable 
millimeter  wave  frequency,  even  in  the  absence  of  an  exter¬ 
nal  biasing  magnet. 

There  has  been  some  progress  for  BaM  films  with  the 
c-axis  normal  to  the  film  plane  (NCA).  Song  et  al 6  made 
such  films  by  pulse  laser  deposition  (PLD)  with  a  near  single 
crystal  of  60  GHz  ferromagnetic  resonance  (FMR)  peak-to- 
peak  derivative  linewidth  (A7/pp)  of  16  Oe.f1  The  remanent 
magnetization  (A/r),  however,  was  quite  small.  Excessive 
stress  limited  film  thicknesses  to  0.85  yum  or  so.  Thicker 
NCA  BaM  films  have  been  made  by  a  hybrid  PLD  liquid 
phase  epitaxy  (LPE)  approach.  *  A  thin  high  quality  PLD 
seed  layer  is  first  deposited.  A  second  layer,  up  to  45  yum  in 
thickness,  is  then  grown  by  LPE  methods.  The  nominal 
60  GHz  A //pp  values  were  in  the  30-75  Oe  range.  The  thick 
film  remanence  ratio  Mrl M s  (RR)-values  were  near  zero,  as 
expected  for  high  quality  thick  films.  NCA  BaM  films  have 
also  been  made  by  screen  printing.9'10  These  films  had  mod¬ 
erate  loss,  with  nominal  60  GHz  A//pp  values  in  the 
300-500  Oe  range,  and  a  high  RR  in  the  range  of  0.96 
range. 

The  new  challenge,  however,  is  to  realize  in-plane  c-axis 
(IPCA)  films  with  low  loss  and  high  RR  values.  The  basic 
advantage  of  IPCA  films  is  in  the  higher  ferromagnetic  reso¬ 
nance  frequency.  For  given  magnetization  and  anisotropy 
values,  the  zero  external  field  FMR  frequency  for  RR=  1 
goes  from  w^=|y|(//4-47rA/i)  for  a  NCA  film  to 
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=  \y\\HA(HA+47rMs)  for  an  IPCA  film.  Here,  y  denotes  the 
gyromagnetic  ratio.  For  the  BaM  parameters  cited  above,  the 
wfmr  / wfm  r  enhancement  factor  is  about  1.45.  This  is  sig¬ 
nificant.  In  principle,  the  realization  of  a  high  RR  for  IPCA 
films  is  also  made  easier,  since  the  in-plane  demagnetizing 
factors  for  thin  films  are  close  to  zero  and  self¬ 
demagnetization  effects  are  minimized. 

Progress  for  IPCA  BaM  films  has  been  limited.  Yoon 
et  a/.11,  2  were  able  to  make  such  films  both  on  a-  and 
m-plane  sapphire  substrates.  High  RR  values  above  0.9  were 
obtained  for  PLD  films  on  a-plane  sapphire  substrates  up  to 
7  yum  or  so  in  thickness,  but  the  linewidths  were  large,  in  the 
range  of  kOe.11  LPE  films  on  m-plane  sapphire  appeared  to 
have  local  regions  with  nominal  60  GHz  A//pp  values  below 
100  Oe,  but  the  overall  linewidth  was  much  larger  and  the 
films  had  relatively  low  RR-values.12  Up  until  now,  no  group 
has  been  able  to  produce  IPCA  BaM  films  that  combine  both 
a  high  remanence  and  low  loss. 

Lattice  mismatch  and  surface  chemistry  provide  the  key 
to  successful  epitaxy.  From  the  narrow  FMR  linewidth  in 
Ref.  6,  it  is  clear  that  good  epitaxy  is  achievable  for  NCA 
BaM  films  on  sapphire  substrates  The  inability  to  produce 
low  linewidth  NCA  films  for  thicknesses  above  0.8  yum  or 
so,  however,  also  points  to  the  presence  of  significant  stress. 
From  Refs.  1 1  and  12,  it  is  clear  that  for  IPCA  films,  both  the 
tf-plane  and  the  ///-plane  sapphire  substrate  growth  options 
present  problems.  For  the  present  work,  a-plane  substrates 
were  used,  but  hybrid  PLD/LPE  films  were  prepared.  Rela¬ 
tively  thin  PLD  Ba  M  epitaxial  films  were  first  grown  on  the 
a-plane  substrate.  This  was  done  in  anticipation  of  a  good 
RR  ratio  for  thicker  LPE  films  grown  on  the  PLD  seed  layer 
in  the  second  step.  Following  Refs.  7  and  8.  it  was  also 
expected  that  the  hybrid  PLD/LPE  approach  would  also 
yield  narrow  linewidths.  This  expectation  was  realized,  but 
some  fine-tuning  in  the  process  parameters  was  needed.  It 
was  necessary,  for  example,  to  use  lower  deposition  tempera¬ 
tures  for  the  initial  PLD  seed  layer  than  normally  used  to 
make  the  highest  quality  film.  It  appears  that  some  degrada¬ 
tion  in  the  seed  layer  is  needed  in  order  to  get  adequate  strain 
relief  for  the  subsequent  thick  LPE  layer. 

The  hybrid  approach  produced  IPCA  oriented  ///-plane 
BaM  films  with  high  remanence  and  moderate  loss.  Films  in 
the  5-7  yum  thickness  range  appear  to  offer  the  best 
remanence/linewidth  compromise  At  7  ^um,  one  has  RR 
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FIG.  I.  (Color  online)  1PCA  BaM  film  XRD  26  scans.  Graphs  (a)  and  (b) 
show  results  for  a  I  y^m  thick  PLD  seed  film  and  the  corresponding  7  fj.m 
thick  hybrid  film,  respectively.  The  BaM,  M,  and  S  miller  indices  labels 
identify  the  important  peaks,  with  S  for  the  sapphire  substrate  and  M  for  the 
characteristic  m  plane  peaks.  The  (h)  inset  shows  a  rocking  curve  for  the  M 
(200)  diffraction  line 

-0.84  and  a  59  GHz.  FMR  A Hpp  value  of  250  Oe.  Thicker 
films  show  a  decrease  in  A//pp  but  also  a  decrease  in  rema- 
nence  At  20  yum,  one  has  RR  —0.33  and  A//pp  — 110  Oe. 
Thinner  films  show  higher  remanence  values  but  degraded 
linewidths,  presumably  due  to  stress  effects. 

The  first  step  was  to  lay  down  the  thin  seed  layer.  A 
nominal  1  yum  thick  epitaxial  PLD  seed  layer  was  grown  on 
u-plane  (1120)  sapphire.  The  energy  fluence  and  repetition 
rate  for  lhe  248  nm  KrF  excimer  laser  was  1.7  J/cnr  and 
25  Hz,  respectively.  The  deposition  time  was  25  min,  the 
substrate  temperature  was  870  °C,  and  the  target  substrate 
distance  was  6.6  cm.  These  conditions  gave  bare  unannealed 
PLD  m-plane  films  with  RR  —0.9  and  a  nominal  60  GHz 
A//pp  of  about  1 .2  kOe,  X-ray  diffraction  (XRD)  data 
showed  that  higher  deposition  temperatures  can  give  some¬ 
what  better  epitaxy  but  no  improvement  in  Iinewidth.  It  is  to 
be  emphasized  that  the  XRD  data  showed  the  growth  of 
m-plane  films.  The  large  linewidths  are  likely  due  to  the  large 
m-plane  BaM  u-plane  sapphire  mismatch  induced  strain.  The 
overall  mismatch  is  complicated  because  of  the  overlay  of 
these  different  classes  of  planes. 

The  second  step  is  to  lay  down  the  thicker  LPE  film  on 
the  seed.  A  mixture  of  BaC03,  Fe203,  and  B203  with  mole 
ratio  of  1.0:0.25:0.571,  as  in  Ref.  8,  was  mixed  thoroughly  at 
1250  °C  for  6  h.  The  homogeneous  melt  was  then  cooled  to 
the  desired  LPE  growth  temperature.  A  growth  temperature 
of  860  °C,  as  measured  with  an  external  thermocouple,  was 
found  to  give  good  epitaxial  films.  Lower  temperatures  gave 
mottled  films.  Higher  temperatures  yielded  no  film  growth 
at  all  As  noted  above,  film  thicknesses  ranged  from 
5  to  20  yum.  Thickness  was  controlled  by  the  immersion 
time  for  the  seed  film  in  the  supersaturated  melt.  The  growth 
rate  was  about  0.2  yum  /min.  The  film  thicknesses  were  de¬ 
termined  from  cross  section  scanning  electron  microscopy 
(SEM). 

Graphs  (a)  and  (b)  in  Fig.  1  show  standard  0!  2 0  diffrac¬ 
tion  scans  for  a  nominal  1  yum  thick  PLD  film  and  a  7  yum 
thick  hybrid  PLD/LPE  film,  respectively.  The  inset  in  (b) 
shows  the  M  (200)  rocking  curve.  The  diffraction  peaks  in 
(a)  come  mostly  from  the  (€00)  /?i-plane  scattering  of  the 
hexagonal  structure,  except  for  a  small  peak  from  the  (114) 
plane.  One  can  also  see  the  a-plane  sapphire  “S”  peak.  In  (b), 


FIG.  2.  (Color  online)  Hysteresis  loops  of  the  magnetic  induction  4ttM  vs 
in-plane  static  field  H.  (a)  and  (b)  show  EA  and  hard  axis  field  loops  and  lhe 
EA  loop  on  an  expanded  field  seale,  respectively,  for  the  7  yum  thick  film  in 
Fig.  I .  (e)  shows  the  expanded  EA  loop  for  lhe  20  fxm  film.  Remanence 
ratios  RRs  and  coercive  fields  (Hc)  indicated  in  (b)  and  (c). 

one  sees  the  same  (€00)  w-plane  reflections  as  in  (a).  Here, 
however,  there  is  no  sapphire  peak  and  the  intensities  of  the 
(€00)  peaks  in  (b)  are  about  three  times  higher  than  those  in 
(a).  Tlie  full  width  at  half  maximum  of  the  XRD  rocktng 
curve  is  0.93°.  While  the  film  is  not  perfectly  oriented,  the 
film  is  reasonably  epitaxial.  All  told,  the  XRD  data  under¬ 
score  the  epitaxial  growth  of  1PCA  BaM  PLD/LPE  hybrid  by 
the  recipe  given  above. 

The  graphs  in  Fig.  2  shows  representative  room  tempera¬ 
ture  vibrating  sample  magnetometry  hysteresis  loop  data  in 
an  induction  4ttM  versus  static  field  H  format.  The  data  in 

(a)  and  (b)  are  for  the  same  7  yum  film  as  in  Fig.  1 .  (a)  shows 
data  with  the  field  H  applied  along  the  in -plane  (IP)  easy  and 
hard  directions,  as  indicated,  (b)  shows  the  easy  axis  (EA)  IP 
loop  on  an  expanded  field  scale,  (c)  shows  similar  data  as  in 

(b)  but  for  a  20  yum  film.  The  data  in  (a)  show  that  the  film 
has  an  in-plane  uniaxial  axis.  The  7  yum  film  saturation  in¬ 
duction  47 tMs  and  effective  uniaxial  anisotropy  field  HA  val¬ 
ues  are  about  3.9  kG  and  17.3  kOe,  respectively.  These  are 
close  to  the  reported  values  for  single  crystal  NCA  BaM 
films.  The  corresponding  AirM  and  //4  values  for  the  20  yum 
thick  him  were  4.0  kG  and  17,3  kOe,  respectively.  The 
20  yum  thick  film  data  in  (c)  shows  much  lower  coercive 
field  and  RR-values  than  do  the  7  yum  film  data  (b). 

Graphs  (a)  and  (b)  in  Fig.  3  show  in-plane  EA  field 
59  GHz  FMR  data  (open  circles)  and  Lorentzian  derivative 
fits  (solid  curves)  for  the  same  films  as  in  Fig.  2,  shown  in  an 
absorption  derivative  versus  field  format.  The  fits  are  reason¬ 
ably  good.  The  FMR  fields  for  the  (a)  and  (b)  data,  taken  as 
the  zero  crossing  points  from  the  fits,  are  about  2.07  and 
2.24  kOe,  respectively.  The  derivative  linewidths  are  as  indi¬ 
cated  in  the  figure  For  IP  EA  field  ferromagnetic  resonance, 
the  applicable  FMR  frequency  (to)  field  ( H0 )  connection  is 
given  by 


= I  y|  v  (//(,  +  Ha)(H  o  +  Ha  +  4  7 tM.)  ,  ( 1 ) 

Based  on  the  operating  frequency  of  59  GH7,  |y|/27r 
=  2.8  GHz/kOe,  and  the  AttM  and  HA  values  cited  above, 
one  obtains  theoretical  FMR  field  values  of  1.91  kOe  (7  yum 
film)  and  1.87  kOe  (20  yum  film).  The  matches  with  the  data 
are  reasonable  and  well  within  the  range  of  error  set  by  the 
uncertainties  in  4itM  and  H A. 
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FIG.  3.  Graphs  (a)  and  (b)  show  59  GHz  FMR  absorption  denvative  vs  H 
profiles  for  the  7  and  20  /xni  thick  films  featured  in  Fig  2,  respectively  The 
open  circles  show  the  data  The  solid  curves  show  Lorentzian  derivative  best 
fits.  The  peak-to-pcak  derivative  lincwidths  (A/Z^,)  are  as  indicated. 


Turn  now  to  the  linewidths.  Recall  that  the  initial  seed 
layer  alone  has  A7/pp~1.2  kOe.  The  250  Oe  linewidth  for 
the  7  (im  hybrid  film  is  about  a  factor  of  five  smaller  and  the 
1 10  Oe  linewidth  for  the  20  /nm  hybrid  film  is  more  than  a 
factor  of  10  lower  than  the  seed  linewidth.  It  is  clear  that  the 
thicker  the  hybrid  film,  the  better  the  crystal  alignment,  the 
more  the  strain  relief,  and  the  better  the  overall  film  quality 
as  measured  by  the  linewidth. 

Figure  4  shows  (a)  RR  and  (b)  59  GHz  A Hpp  linewidth 
values  as  a  function  of  film  thickness.  From  (a),  one  sees  that 
the  RR  value  decreases  more  or  less  smoothly  from  about  0.9 
at  1  fim  down  to  about  0.3  at  20  /im.  As  discussed  in  con¬ 
nection  with  Fig.  2,  the  decrease  in  the  RR  accompanied  by 
a  decrease  in  the  coercive  force  as  well.  These  decreases  are 
tentatively  ascribed  mainly  to  the  strain  relief  that  comes 
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FIG  4  (a)  RR  vs  film  thickness,  (b)  59  GHz  peak-to-peak  derivative  line- 
width  A/Zpp  vs  film  thickness.  The  data  points  for  the  l  film  are  a  pure 
PLD  film,  with  the  FMR  measurements  done  at  about  66  GHz. 


with  thicker  LPE  layers  and  the  corresponding  reduction  in 
local  stress  related  obstacles  to  domain  wall  movement. 
Graph  (b)  shows  a  similar  trend  for  the  linewidth.  The  nomi¬ 
nal  Vband  linewidth  decreases  from  about  1200  Oe  at  1  ^m 
down  to  about  1 10  Oe  for  20  (im  thickness.  Roughly  speak¬ 
ing,  one  can  see  that  a  thickness  around  7  /zm  or  so  offers  a 
reasonable  compromise  for  a  high  RR  and  low  linewidth, 
with  RR  —0.84  and  A/7pp  — 250  Oe.  Thinner  films  offer  high 
RR  and  high  loss,  while  thicker  films  offer  small  RR  and  low 
loss. 

In  summary,  hybrid  PLD  and  LPE  techniques  have  been 
used  to  produce  1PCA  BaM  films  with  low  microwave  loss 
and  a  high  remanence.  The  use  of  an  initial  thin  PLD  seed 
film  of  moderate  quality  with  a  LPE  grown  thicker  film  on 
top,  appears  to  provide  a  viable  way  to  produce  useful  films 
for  planar  high  frequency  MM1C  devices.  It  is  noteworthy, 
moreover,  that  the  llOOe  derivative  linewidth  for  the 
20  /im  film  is  one  of  the  smallest  values  achieved  lor  IPCA 
oriented  barium  ferrite  film  on  any  substrate. 

This  work  was  supported  in  part  by  the  Oftice  of  Naval 
Research  (ONR)  Grant  Nos.  N0001 4-07- 1-0597  and 
N000 14-08- 1-1 050  (through  sub-award  No.  PT 103701- 
SC 10 1157  from  Virginia  Commonwealth  University),  the 
ARO-DARPA  Multifunctional  Materials  Seedling  Program, 
Contract  No.  W91  INF-06-1-0163,  and  ARO-MUR1  Grant 
No.  W911NF-04- 1-0247  Dr.  S.  Kohli  and  Dr.  P  McCurdy 
are  acknowledged  for  the  XRD  and  SEM  measurements. 

‘S.  Capraro,  T.  Rouiller.  M.  L.  Berre,  J.  P.  Chatelon,  B.  Bayard,  O.  Barbier. 
and  J.  J.  Rousseau.  IEEE  Trans.  Compon  Packag.  Technol.  30.  411 
(2007). 

2C.  K.  Qucck  and  L.  E.  Davis,  Electron.  Lott.  39,  1595  (2003) 

3P.  Kwan  and  C.  Vittoria,  IEEE  Trans.  Magn.  29,  3431  (1993). 

4P.  Kwan  and  C.  Vittoria.  IEEE  Trans  Microwave  Theory  Tech.  41,  652 
(1993). 

5C.  A.  Carosella,  D  B.  Chrisey,  P  Lubitz,  J  S.  Horwitz,  P  Dorsey,  R. 
Seed,  and  C.  Vittoria,  J.  Appl.  Phys.  71,  5107  (1992). 

6Y.  Y  Song,  S  Kalarickal,  and  C.  E  Patton.  J  Appl  Phys  94,  5103 
(2003). 

7S.  D.  Yoon  and  C.  Vittoria,  J.  Appl.  Phys.  96,  2131  (2004). 

*S  G  Wang,  S.  D.  Yoon,  and  C.  Vittoria,  J  Appl.  Phys.  92,  6728  (2002). 
9V.  G.  Hams,  Z.  Chen,  Y.  Chen,  S.  Yoon,  T  Sakai,  A.  Gielcr,  A  Yang.  Y. 
He,  K.  S.  Zicmcr.  N.  X.  Sun,  and  C.  Vittona,  J  Appl.  Phys.  99,  08M91I 
(2006). 

I0Y.  Chen,  T  Sakai,  T.  Chen,  S.  D  Yoon,  A.  L.  Geiler.  C.  Vittoria.  and  V.  G. 

Hams,  Appl  Phys.  Lett  88  062516  (2006). 
nS.  D.  Yoon,  C.  Vittoria,  and  S  A.  Oliver,  J  Appl.  Phys.  93,  4023  (2003). 
,2S  D.  Yoon  and  C.  Vittona,  J.  Appl.  Phys.  93,  8597  (2003). 


Author  complimentary  copy.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://apl.aip.org/apl/copyrightjsp 


Das,  el  aif  ‘'Electric  field  tunable  ...  heterostructures  ...”  Distribution  version 


ADVANCED 

/MATERIALS 


Electric  Field  Tunable  Low  Loss  Multiferroic  Ferrimagnetic-ferroelectric  Heterostructures[*] ** 

By  Jaydip  Das ,  Young- Yeal  Song,  Nan  Mo,  Pavol  Krivosik,  and  Carl  E.  Patton 
(Distribution  version,  accepted  for  publication  in  Advanced  Materials,  12  December  2008) 


Multiferroic  materials  have  attracted  major  attention  in  recent 
years  due  to  the  unique  possibility  to  tune  the  magnetic 
properties  with  a  modest  electric  field  and  vice  versa.11'31  Up 
to  now,  however,  no  single  phase  material  has  been  put 
forward  that  demonstrates  a  practical  capacity  for  such  tuning 
at  room  temperature.141  This,  in  turn,  has  resulted  in  a  shift  in 
focus  from  single  phase  multiferroics  to  multiferroic 
heterostructures. [<1  Besides  the  new  materials  aspect  of  such 
structures,  there  are  numerous  device  benefits.  These  include 
fast  wide  bandwidth  tuning  with  low  insertion  loss  in  the 
microwave  and  millimeter  wave  frequency  range,  planar 
integratibility,  and  negligible  power  consumption.  These 
benefits  speak  to  the  strong  potential  for  practical  information 
storage,  logic  device,  radar,  satellite,  and  telecommunication 
applications. 

Both  multilayer  thin  film  stacks  as  well  as  composite 
nanopillar  films  can  be  produced  by  various  methods  to  satisfy 
such  needs.16  81  The  layer-by-layer  fabrication  of  multilayer 
structures  is  generally  more  straightforward  than  the 
complicated  growth  process  needed  for  nanopillar  films. 
Previous  workers  have  studied  a  variety  of  layered  stacks  with 
a  broad  range  of  ferromagnetic  or  ferrimagnetic  (FM)  and 
ferroelectric  (FE)/piezoelectric  (PE)  components. {9  l3]  Up  to 
now,  the  main  emphasis  has  been  on  FM-PE  layers  and  the 
tuning  of  the  magnetic  response  through  the  electric  field 
induced  stress  on  the  FM  layer.  In  our  view,  however,  electric 
field  tunable  hybrid  magneto-electric  modes  in  the  monolithic 
FM-FE  structures  can  yield  a  better  tuning  of  the  magnetic 
response  with  electric  fields  and  a  better  compatibility  with 
various  device  applications. 

In  principle,  a  micron  thick  FM-FE  monolithic 
heterostructure,  in  combination  with  10  nm  or  so  embedded 
biasing  electrodes,  for  example,  can  offer  (1)  a  good  quasi- 
lattice-to-lattice  contact  with  the  layers  to  give  rise  to  the  best 
possible  magneto-electric  coupling  and  (2)  an  appreciable 
magnetic  tuning  at  small  applied  voltages  in  the  range  of  15-25 
V  or  so.  However,  the  challenge  to  fabricate  heteroepitaxial 
ferrimagnetic  oxide  -  metal  electrode  -  ferroelectric  oxide 
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stacks  with  the  key  physical  properties  of  each  layer  properly 
maintained  and  with  extremely  smooth  interfaces,  has 
restricted  the  realization  of  such  stacks  so  far.  There  have 
been  some  limited  works  with  thick  mechanically  assembled 
structures  that  (1)  does  not  offer  the  lattice-to- lattice  contact 
between  the  layers  and  (2)  needs  about  two  orders  of 
magnitude  higher  voltages  than  the  monolithic  structures.114, 151 

In  this  report,  we  present  the  realization  of  a  significant 
tuning  of  magnetic  response  with  extremely  small  applied 
voltages  in  a  multiferroic  FM-FE  heterostructure.  The 
structure  was  comprised  of  a  pulsed  laser  deposited  near  single 
crystal  yttrium  iron  garnet  (YIG)  layer  with  very  low 
microwave  loss,  an  oriented  barium  strontium  titanate  (BSTO) 
layer  with  good  electric  field  tunable  relative  dielectric 
permittivity  en  and  embedded  platinum  electrodes  on  a  single 
crystal  gadolinium  gallium  garnet  (GGG)  substrate.  The  results 
indicate  an  order  of  magnitude  higher  magneto-electric 
coupling  than  for  pasted  non-monolithic  structures. 

Scheme  1  shows  a  concept  diagram  of  the  cube-on-cube 
unit  cell  match-up  for  an  all  oriented  YIG-Pt-BSTO-Pt-GGG 
stack.  The  cubic  YIG  lattice,  with  a  lattice  parameter  of  12.38 
A,  matches  to  the  GGG  lattice  within  0  03%.  This  close 
match,  along  with  the  compatible  oxygen  coordinations,  serves 
to  promote  epitaxial  YIG  film  growth  The  lattice  parameter 
for  platinum  is  about  3.93  A.  As  apparent  from  the  diagram, 
each  YIG  unit  cell  can  reasonably  accommodate  nine  face- 
centered-cubic  (FCC)  Pt  unit  cells  with  a  mismatch  of  about 
5%.  Moving  up  the  stack,  the  FCC  BSTO  lattice  parameter  is 
3.95  A,  for  a  mismatch  of  about  0.5%  with  the  Pt  All  told, 
this  overall  assembly  is  expected  to  yield  a  reasonable 
heteroepitaxial  layered  structure  For  ease  of  visualization,  the 
concept  diagram  shown  here  is  for  (100)  plane  layers.  Similar 
matchups  apply  to  the  actual  structure  developed  here  with 
(111)  plane  heteroepitaxy. 

The  initial  YIG  layer  was  deposited  in  a  100  mTorr 
oxygen  atmosphere  on  the  (111)  GGG  substrate  heated  to 


Scheme  1.  Concept  diagram  of  the  cube-on-cube  unit  cells  of  the 
YIG-Pt-BSTO-Pt  heterostructure  on  the  GGG  substrate.  Each  box 
represents  one  unit  cell  of  the  respective  material,  as  indicated.  The 
atoms  in  the  respective  unit  cells  arc  also  shown. 
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about  800  °C.  The  thickness  was  held  at  a  nominal  value  of 
about  0.3  pm.  This  is  about  as  thick  as  one  can  go  and  still 
avoid  a  degradation  in  film  quality  and  increase  in  microwave 
loss  due  to  the  lattice  mismatch  induced  strain.  This  phase  of 
the  structure  development  is  relatively  straightforward.  The 
challenge  is  to  achieve  near  single  crystal  films  with  narrow 
ferromagnetic  resonance  (FMR)  linewidths.  The  bare  film, 
without  the  processing  to  obtain  the  additional  layers,  showed 
a  9.5  GHz  derivative  FMR  linewidth  of  about  2  Oe.  The 
processing  needed  to  produce  the  subsequent  layers  provides 
some  annealing  of  the  YIG  layer  and  a  linewidth  reduction  to 
about  0.7  Oe. 

The  most  critical  step  in  the  heteroepitaxial  multilayer 
development  lies  in  the  BSTO  and  embedded  Pt  electrode 
layer  growth  on  top  of  the  YIG  layer.  The  first  step  is  to  lay 
down  a  thin,  epitaxial,  continuous,  and  pinhole  free  Pt  layer  to 
serve  as  the  seed  for  the  oriented  and  continuous  BSTO  layer 
that  comes  next.  The  Pt  layer  must  be  thin,  on  the  order  of 
tens  of  nm,  so  that  the  hybrid  YIG-BSTO  mode  tuning  is  not 
unduly  compromised  by  the  electrode  in  between  In  order  to 
realize  an  epitaxial  Pt  layer,  high  temperature  processing  is 
needed  at,  for  example,  600  °C  or  so.  However,  the  difference 
in  the  diffusion  kinetics  between  oxides  (YIG)  and  metals  (the 
Pt)  at  such  temperatures  generally  leads  to  an  island  like 
morphology  that  leaves  the  oxide  surface  partially  exposed.*1 
This  is  especially  problematic  for  the  requisite  thin  Pt  layers 
needed  in  this  work.  A  special  strategy  was  used  here.  The 
first  step  was  to  deposit  a  continuous  but  polycrystalline  Pt 
layer  with  the  GGG-YIG  base  structure/substrate  held  at  room 
temperature.  The  subsequent  high  temperature  BSTO 


deposition  step,  considered  in  more  detail  below,  provide  an 
anneal  for  the  platinum  layer.  This  anneal  produced  the 
desired  (111)  oriented  state.  Note,  however,  that  such  a  post 
deposition  anneal  can  also  lead  to  island  formation.  For  this 
reason,  it  was  necessary  to  lay  down  a  fairly  thick  30  nm  layer 
of  Pt  in  order  to  ensure  a  continuous,  epitaxial,  (111)  oriented 
layer.  This  thickness  is  still  about  two  orders  of  magnitude 
smaller  than  the  platinum  skin  depth  at  10  GHz,  for  example, 
so  that  the  desired  magneto-electric  coupling  between  the  YIG 
and  the  BSTO  is  not  substantially  affected. 

Consider  next  the  ferroelectric  BSTO  layer  and  the  top  Pt 
electrode  capping  layer.  Two  structures  were  prepared  with 
nominal  1  and  5  pm  thick  BSTO  layers.  These  layers  were  put 
down  at  650  °C  in  75  mTorr  of  oxygen.  There  was  one  caveat. 
Direct  BSTO  deposition  on  the  oriented  Pt  gives  rise  to  a 
polycrystalline  layer  with  numerous  cracks.  In  order  to  obtain 
a  (111)  crack  free  BSTO  layer,  a  very  thin  nominal  5-10  nm 
MgO  buffer  layer  was  first  deposited  on  the  Pt.  The  1  and  5 
pm  BSTO  thickness  choices  were  based  on  several  trade-offs. 
(1)  The  thinner  the  layer,  the  less  the  applied  voltage  needed  to 
get  the  same  electric  field  across  the  FE  layer.  (2)  Thin  layers, 
however,  also  mean  a  reduced  magneto-electric  hybrid  mode 
coupling  between  the  FM  and  FE  components.  For  the 
selected  thicknesses,  a  20  V  BSTO  layer  tuning  voltage 
corresponds  to  an  electric  field  in  the  106  V/m  range.  This 
should  be  sufficient  to  give  strong  magneto-electric  coupling 
and  reasonable  tuning  of  the  magnetic  response. 141  Finally,  a 
50  nm  Pt  layer  was  deposited  on  BSTO  at  room  temperature  to 
complete  the  YIG-Pt-BSTO-Pt  structure 

X-ray  diffraction  (XRD),  along  with  combined  scanning 
electron  microscopy  (SEM)  and  energy  dispersive 
spectroscopy  (EDS)  measurements,  were  used  for  the 
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Figure  1.  (a)  X-ray  diffraction  20  scan  for  the  heterostructure,  (b)  Element  specific  energy  dispersive  spectrum  (EDS)  results  for  the  0.3  pm 
thick  YIG  and  the  1  pm  thick  BSTO  layered  structure  with  embedded  Pt  electrodes,  as  indicated.  The  data  have  been  taken  from  the  side,  (c) 
Scanning  electron  micrograph  (SEM)  of  the  BSTO  surface.  The  purple  triangle  marks  one  (100)  facet,  (d)  Schematic  BSTO  unit  cell  diagram 
with  the  (1 1 1)  plane  (marked  as  grey)  oriented  horizontally,  with  one  (100)  facet  marked  with  a  purple  triangle,  as  in  (c). 
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structural,  surface,  and  interface  characterizations.  Hysteresis 
loop  and  FMR  measurements  were  used  for  the  magnetic 
characterizations.  Capacitance  measurements  were  used  to 
determine  the  relative  permittivity  £r  and  the  voltage  tunability 
of  the  BSTO.  Finally,  the  FMR  response  and  the  tuning  of  this 
response  were  measured  as  a  function  of  the  applied  BSTO 
layer  bias  voltage.  The  data  indicate  good  (111) 
crystallographic  orientation  for  all  layers,  a  low  FMR 
linewidth  for  the  YIG,  a  high  er  and  good  voltage  tunablility 
for  the  BSTO,  and  a  strong  magneto-electric  hybrid  mode 
coupling  between  the  YIG  and  the  BSTO  layers. 

Figure  1  shows  representative  XRD,  EDS,  and  SEM 
results.  Graph  (a)  shows  the  XRD  pattern  and  schematic 
stacking  diagram  for  the  heterostructure.  The  origin  of  the 
various  peaks  are  as  labeled.  One  sees  an  intense  YIG  444 
peak  that  establishes  the  nearly  ideal  (1 11)  orientation  of  the 
layer.  The  diagram  also  show  weaker  but  distinct  peaks 
associated  with  the  epitaxial  BSTO  (111)  and  Pt  (111)  layers. 
The  MgO  is  not  discernable  from  the  XRD  data  in  (a)  due  to 
the  extreme  thinness  (5-10  nm)  of  the  layer.  Thicker  MgO 
layers  (e.g.,  100  nm)  do  show  clear  (111)  XRD  peaks. 

Panel  (b)  shows  the  representative  EDS  mapping  for  the 
structure  with  the  1  pm  thick  BSTO  layer,  in  a  side  view 
perspective.  The  thicknesses  of  visible  YIG,  Pt,  and  BSTO 
layers  are  as  labeled.  A  composition  analysis  of  the  YIG  and 


BSTO  layers  confirms  the  expected  atom  fractions  for  both  the 
materials.  Panel  (c)  shows  a  SEM  micrograph  of  the  BSTO 
surface  One  can  see  pyramidal  shaped  BSTO  grains  that  point 
out  of  plane.  For  a  given  grain,  one  side  of  the  pyramid 
corresponds  to  a  (100)  facet.  One  such  facet  is  marked  by  a 
purple  triangle  in  (c).  Diagram  (d)  shows  the  schematic  layout 
of  the  BSTO  unit  cell.  One  of  the  horizontal  (111)  planes  is 
marked  in  grey.  From  this  diagram,  one  can  see  clearly  the 
pyramidal  shape  evident  from  the  micrograph  in  (c)  and  the 
associated  (100)  facets.  The  SEM  results,  therefore,  serve  to 
confirm  directly  the  (1 1 1)  BSTO  orientation  given  by  the  XRD 
data. 

As  mentioned  earlier,  the  very  good  match  between  the 
YIG  and  GGG  lattice  is  mainly  responsible  for  the  nearly 
single  crystalline  growth  of  the  YIG  layer  along  the  (111) 
direction.  The  high  melting  point  and  high  surface  energy 
anisotropy  then  help  to  orient  the  embedded  Pt  layer  along  the 
most  favorable  (111)  direction  during  the  high  temperature 
anneal  that  comes  from  the  subsequent  MgO  and  BSTO 
depositions. [1  1  The  (111)  MgO  buffer  serves  to  promote  the 
growth  of  (111)  BSTO.  Note  that  the  BSTO  layer,  when 
deposited  directly  on  the  Pt  layer,  is  polycrystalline  in  nature 
and  shows  cracks.  This  indicates  a  reasonable  amount  of  strain 
on  the  BSTO  from  the  5%  lattice  mismatched  YIG-Pt 
interface.  Most  likely,  the  6%  bigger  MgO  unit  cell  relative  to 


Figure  2.  Magnetic  and  electric  response  results  for  the  heterostructure.  Hysteresis  loop  data  for  (a)  in-plane  (IP)  and  (b)  perpendicular-to- 
plane  (PTP)  field  orientations.  Ferromagnetic  resonance  absorption  derivative  vs.  static  field  H  spectra  at  9.5  GHz  for  (c)  IP  and  (d)  PTP 
field  orientations.  Labels  designate  the  main  mode  and  the  higher  order  spin  wave  modes,  (e)  Dielectric  constant  Er  vs.  electric  field  E 
applied  across  the  BSTO  layer.  The  closed  and  open  circles  show  the  data  for  positive  to  negative  and  negative  to  positive  field  sweeps,  as 
indicated. 
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the  platinum,  serves  to  compensate  the  YIG  induced  strain  and 
allow  for  an  oriented  and  crack-free  BSTO  layer.  More 
extensive  crystal  chemistry  related  experiments  coupled  with 
theoretical  analyses  will  be  needed  to  obtain  a  better  insight  to 
the  complicated  nature  of  these  oxide-metal-oxide  interfaces. 

Figure  2  shows  representative  room  temperature  hysteresis 
loop,  FMR,  and  dielectric  response  results  for  the  oriented 
heterostructure.  Graphs  (a)  through  (d)  characterize  the 
magnetic  properties  of  the  YIG  layer.  Graphs  (a)  and  (b)  show 
hysteresis  loops  in  a  magnetic  induction  4 n\i  vs.  static  field  H 
format  for  in-plane  (IP)  and  perpendicular-to-plane  (PTP) 
static  fields,  respectively.  Note  the  different  field  scales  for 
the  IP  and  PTP  loop  data.  Graphs  (c)  and  (d)  show  the  9.5 
GHz  IP  and  PTP  field  FMR  responses,  respectively,  in  a 
normalized  absorption  derivative  vs.  H  format.  The  strongest 
peaks  are  labeled  as  the  main  mode  in  both  graphs.  The  less 
intense  peaks  on  the  low  field  side  of  the  spectra  in  graph  (d) 
are  labeled  as  higher  order  spin  wave  modes.  Graph  (e)  shows 
the  BSTO  relative  dielectric  constant  sr  vs.  the  dc  electric  field 
E  applied  across  the  layer  as  obtained  from  the  100  kHz 
capacitance  measurements. 

Consider  first  the  hysteresis  loop  data  in  (a)  and  (b).  For 
the  IP  field  case  in  (a),  one  sees  a  nearly  square  loop  with  a 
very  low  coercive  force.  The  nominal  saturation  induction 
4 7rMs  of  about  1.7  kG  is  close  to  literature  values  for  single 
crystal  YIG.1181  For  the  PTP  field  case  in  (b),  the  loop  is 
definitely  not  square  and  saturation  is  realized  only  for  fields 
above  1800  Oe  or  so,  close  to  the  saturation  field  expected  for 
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Figure  3.  Electric  field  tuning  of  the  magnetic  response  for  the 
he tero structure  (a)  The  FMR  absorption  derivative  vs.  in-plane 
static  field  H  spectra  at  9.5  GHz  for  0  and  25  V  applied  voltages, 
as  indicated.  Inset:  The  FMR  resonance  position  as  a  function  of 
the  applied  voltage  V.  (b)  Incremental  frequency  shift  vs.  applied 
voltage.  The  circles  show  the  data  and  the  solid  line  is  a  guide  to 
the  eye. 


an  isotropic  YIG  film.  This  indicates  that  the  YIG 
magnetization  lies  in  the  film  plane,  a  state  that  is  expected  to 
yield  the  strongest  magneto-electric  coupling  for  the  Y1G- 
BSTO  heterostructure. 

Consider  now  the  FMR  results  in  (c)  and  (d).  The 
relatively  low  and  high  main  mode  fields  for  the  IP  and  PTP 
field  configurations  are  consistent  with  the  hysteresis  loop 
properties  in  (a)  and  (b).  The  very  narrow  0,7  Oe  main  mode 
derivative  linewidth  for  the  PTP  field  case  in  (d)  is  only 
slightly  larger  than  the  0.5  Oe  or  so  near-intrinsic  linewidth 
found  for  bare  YIG  films.  The  presence  of  weak  higher  order 
spin  wave  modes  at  fields  below  the  main  mode  gives 
additional  evidence  of  the  layer  quality.  The  somewhat  larger 
3  Oe  linewidth  for  the  IP  field  case  in  (c)  is  likely  due  to  a 
small  contribution  from  the  inhomogeneity  related  two 
magnon  scattering  in  the  YIG  layer. 

Finally,  consider  the  electrical  response  properties  of  the 
BSTO  layer  in  (e).  Overall,  one  sees  a  significant  tuning 
response  for  the  relative  dielectric  permittivity  er  with  the 
electric  field  E.  The  fact  that  the  response  shows  some  small 
amount  of  hysteresis  points  to  the  ferroelectric  nature  of  the 
BSTO  layer.  The  cusp  response  and  the  nearly  35%  change  in 
er  at  15  V/pm  demonstrate  the  very  good  electric  field 
tunability  of  the  layer.  The  value  of  er  at  E  »  0  is  close  to 
literature  values  for  the  BSTO  films  1191  The  important 
message  from  these  data  is  that  (1)  the  magnetic  responses  are 
consistent  with  the  known  properties  of  near  perfect  YIG 
single  crystals,  (2)  the  electric  response  shows  the  presence  of 
a  viable,  high  quality,  high  en  and  tunable  BSTO  layer,  and  (3) 
as  such,  the  high  quality  YIG  and  BSTO  layers  realized  here 
are  expected  to  be  amenable  to  a  strong  hybrid  mode  electric 
field  tuning  response. 

Figure  3  shows  the  change  in  the  magnetic  response  due  to 
the  application  of  electric  fields  for  the  heterostructure  with  the 
5  pm  BSTO  layer.  The  main  graph  in  (a)  shows  the  9.5  GHz 
IP  field  FMR  profile  in  an  absorption  derivative  vs.  static 
magnetic  field  (If)  format  for  applied  voltages  of  zero  and  25 
V  across  the  BSTO  layer,  as  indicated.  The  inset  shows  the 
FMR  resonance  position  vs.  applied  voltage  V  profile.  Both 
sets  of  data  are  for  increasing  voltage.  Graph  (b)  shows  data  in 

(a) ,  but  shown  in  terms  of  the  incremental  frequency  shift  in 
MHz/V  as  a  function  of  voltage.  The  structure  with  the  1  pm 
BSTO  layer  showed  no  observable  electric  field  driven  FMR 
shift. 

From  graph  (a),  one  sees  that  25  V  causes  the  FMR 
response  to  shift  up  in  field  by  about  5  Oe.  From  the  inset 
data,  one  sees  a  small  down-shift  to  a  shallow  minimum  at 
about  5  V  followed  by  an  incremental  shift  that  increases 
rapidly  with  V.  The  incremental  frequency  response  in  graph 

(b)  follows  the  same  trend  as  in  the  inset  of  (a).  The  overall 
shape  of  this  response  is  as  one  would  expect  for  the  hysteretic 
ferroelectric  BSTO  material. 

The  electric  field  tuning  of  the  magnetic  resonance 
revealed  by  the  data  in  Fig.  3  is  very  substantial.  Taken  at  face 
value,  the  5  Oe  FMR  field  shift  in  (a)  amounts  to  an  equivalent 
frequency  shift  of  14  MHz  and  the  corresponding  tunability  is 
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about  0.6  MHz/V.  This  is  about  two  times  the  hybrid  mode 
tunability  reported  in  Ref.  [15]  for  mechanically  pasted  non- 
monolithic  structures  at  a  voltage  of  400  V.  The  incremental 
frequency  response  in  (b)  shows  an  even  higher  tunability'.  As 
evident  from  the  graph,  the  incremental  frequency  tuning  is 
about  2  MHz/V  at  25  V  This  is  nearly  one  order  of  magnitude 
greater  than  that  reported  previously.111  The  good  quasi- 
lattice-to-lattice  contact  is  likely  responsible  for  this  enhanced 
tunability.  This,  in  turn,  is  indicative  of  the  very  good 
magneto-electric  coupling  between  the  layers.  Moreover,  the 
20  -  25  V  bias  voltages  needed  for  the  good  coupling  can  be 
achieved  easily  in  such  multilayered  structures. 

There  are  two  further  points  for  consideration.  First,  for 
these  monolithic  structures  as  fabricated,  at  least  some  of  the 
response  in  Fig.  3  may  be  due  to  stress  tuning.  Second,  in  the 
absence  of  a  full  theoretical  analysis  of  the  hybrid  mode 
response,  it  is  not  yet  clear  whether  the  combined  hybrid  mode 
and  stress  tuning  effects  should  add  or  subtract.  A  full 
theoretical  understanding  of  these  processes  should  allow  for 
the  full  optimization  of  the  tuning  and  realize  an  even  higher 
tuning  response  than  reported  here  These  considerations  make 
the  monolithic  heterostructure  approach  extremely  promising 
for  electric  field  controlled  small  and  fast  magnetic  devices. 
The  ability  to  accomplish  such  control  in  planar  monolithic 
structures,  as  demonstrated  here,  holds  great  promise  for  many 
practical  applications  in  RF  electronics,  communications, 
electronic  logic,  and  radar  technology. 

In  summary,  an  electric  field  tunable  magneto-electric 
monolithic  multilayered  heterostructure  with  ferrimagnetic  and 
ferroelectric  oxide  layers,  along  with  embedded  metal 
electrodes,  has  been  achieved.  The  “all  oriented” 
heteroepitaxial  stack  shows  the  desirable  physical  properties  of 
a  ferrimagnetic  layer  with  very  low  microwave  loss  and  a 
ferroelectric  layer  with  high  permittivity  and  good  electric 
field  tunability.  The  structure  shows  a  significant  electric  field 
tunability  of  the  coupled  hybrid  magneto-electric  modes  for 
small  applied  voltages.  This  response  is  an  order  of  magnitude 
better  than  previously  obtained  results  for  thick  mechanically 
bonded  layers. 

The  realization  of  such  a  hybridized  structure  is  the  first 
critical  step  in  the  fabrication  of  practical  voltage  tunable 
magnetic  devices.  Future  work  should  focus  on  further 
increases  in  the  tuning  response  through  structure  optimization 
and  the  use  of  multiple  magnetic  and  electric  layers.  In  the 
broad  view,  the  materials  advances  realized  here  for  the  high 
temperature  growth  of  high  quality  ferromagnetic  oxide-metal- 
ferroelectric  oxide  stacks  with  very  strong  magneto-electric 
coupling  is  expected  to  have  a  significant  impact  on  the 
physics,  chemistry,  and  materials  science  of  the  thin  film 
heterostructures. 
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The  frequency  dependence  of  the  ferromagnetic  resonance  linewidth  AH  in  polycry  stalline, 
hot  isostatically  pressed  pure  and  Ca-Ge  substituted  yttrium  iron  garnet  (YIG)  spheres  was 
measured  between  1.95-6  GHz  and  analyzed.  The  linewidths  for  these  ultradense  ferrite  materials 
arein  the  5-40  Oe  and  5-25  Oe  ranges  for  the  pure  and  the  substituted  YIG  samples  respectively. 
The  A H  vs.  frequency  data  show  especially  sharp  Buffler  peaks  due  to  the  band  edge  crossover 
effect  that  is  the  signature  for  grain-to-grain  (GG)  two  magnon  scattering  (TMS)  [J.  Appl  Phys  30, 
172S  (1959)].  Both  the  GG  and  recent  grain  boundary  (GB)  TMS  models  were  used  to  obtain 
quantitative  fits  to  the  data.  The  fits  demonstrate  the  expected  GG  TMS  contribution  to  the 
linewidth  for  frequencies  above  the  Buffler  peak.  They  also  reveal  two  unexpected  results,  (1)  a 
dominant  GB  TMS  role  for  frequencies  below  the  Buffler  peak  where  the  GG  TMS  linewidth 
contribution  is  essentially  zero,  and  (2)  a  significant  GB  TMS  role  for  frequencies  above  the  peak 
where  the  linewidth  has  generally  been  associated  with  GG  processes  only.  In  the  high  frequency 
regime,  the  GB  scattering  term  appears  to  explain  about  half  of  the  overall  linewidth. 

PACS  number(s):  76.50 +g,  75.30.Ds,  75.50.Gg, 


I.  INTRODUCTION 

Despite  many  decades  of  study,  the  origins  of 
microwave  loss  in  polycrystalline  magnetic  materials 
remain  a  hot  topic  of  research.  There  have  been  recent 
advances  on  several  fronts,  made  possible  in  a  large  part  by 
the  availability  of  ultra  dense  ferrite  samples  with  near  zero 
porosity  and  second  phase.1 2  The  notable  work  of  Nazarov 
and  co-workers/  for  example,  was  the  first  to  demonstrate 
that  ultra  dense  polycrystalline  ferrites  can  show  close-to- 
theoretical  linewidths  based  on  grain-to-grain  (GG)  two 
magnon  scattering  (TMS). 

More  recent  work  of  Mo  and  co-workers,4  in  addition, 
has  advanced  the  understanding  of  TMS  relaxation 
processes  in  two  ways.  First,  the  application  of  clever 
metrology  methods  led  to  order  of  magnitude  improvements 
in  the  accuracy  of  high  field  effective  linewidth 
measurements.  This  has  led,  in  turn,  to  the  discovery  of  a 
new  role  for  hybrid  electromagnetic-exchange  spin  waves 
and  new  grain  boundary'  (GB)  two  magnon  scattering  in 
polycrystalline  ferrites.  Second,  and  in  a  completely 
different  context,  even  more  recent  ferromagnetic  resonance 
(FMR)  linewidth  measurements  in  nano-crystalline  metallic 
magnetic  films  have  demonstrated  the  crucial  roles  of  both 


grain-to-grain  and  grain  boundary  two  magnon  scattering  in 
microwave  relaxation. 

This  work  returns  to  the  problem  of  microwave 
relaxation  in  ultra  dense  polycrystalline  ferrites  The  results 
from  Ref.  3  have  been  extended  to  the  low  GHz  regime  and 
a  comprehensive  analysis  of  the  full  range  of  linewidth  data 
over  the  characteristic  Buffler  peak  in  a  linewidth  vs. 
frequency  data  display  format.  This  peak  occurs  for  sphere 
samples  when  the  FMR  point  moves  above  the  top  edge  of 
the  spin  wave  band  at  low  wave  numbers.  For  the  yttrium 
iron  garnet  (YIG)  materials  used  here,  this  effect  occurs  at 
about  3  GHz.  The  analysis  shows  that  both  GG  and  GB 
processes  play  important  roles  in  the  structure  of  this 
Buffler  peak. 

The  paper  is  organized  as  follows.  Section  II  gives  a 
brief  qualitative  description  of  two  magnon  scattering 
processes  in  bulk  ferrites  and  the  qualitative  origins  of  the 
Buffler  peak.  Section  III  elaborates  sample  properties, 
gives  representative  FMR  profiles,  and  summarizes  the 
overall  linewidth  vs.  frequency  results.  Qualitative 
connections  to  two  magnon  scattering  processes  are  also 
covered  in  this  section.  Section  IV  provides  working 
equations  for  the  quantitative  analysis  of  the  data.  Section 
V  presents  theoretical  fits  to  the  data  and  makes  specific 
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connections  to  the  linewidth  contributions  due  to  different 
scattering  processes.  Section  VI  gives  a  summary  and 
conclusion. 

II.  A  VIGNETTE  ON  TWO  MAGNON 
SCATTERING  IN  BULK  FERRITES 

The  early  FMR  data  for  ferrite  materials  generally 
showed  larger  linewidths  than  expected  for  intrinsic 
processes.  Theoretical  considerations  by  Clogston  et  a/.8 
and  Schlomann  1  showed  two  possible  origins  of  these 
large  linewidths,  TMS  processes  and  simple 
inhomogeneous  linebroadening.  Inhomogeneity  effects 
dominate  when  the  corresponding  local  fields  are  large 
compared  to  the  saturation  induction  AnMs .  Two  magnon 
scattering,  while  still  based  on  inhomogeneities,  is  more 
subtle.  TMS  processes  rely  on  small  local  effective  field 
spatial  variations  due  to  randomly  oriented  crystalline  grain 
in  anisotropic  materials,  dipole  fields  due  to  grain 
boundaries,  second  phase,  pores,  surface  pits, 
magnetostriction,  and  so  on.  These  field  fluctuations  lead 
to  a  coupling  between  the  driven  mode  and  available  spin 
wave  modes  at  the  same  frequency.  The  coupling  provides 
an  additional  channel  for  energy  relaxation  from  the  driven 
mode,  first  to  the  degenerate  spin  waves,  and  eventually  to 
the  thermal  bath.  This  is  the  essence  of  two  magnon 
scattering. 

The  first  experimental  signature  of  two  magnon 
scattering  as  a  valid  FMR  relaxation  process  came  from  the 
seminal  work  of  Buffler  1 959. 7  Schldmann’s  work,  noted 
above,  had  pointed  out  a  peculiar  spin  wave  band  effect  for 
sphere  shaped  samples.  As  one  increases  the  applied  field, 
there  is  a  point  at  which  the  FMR  frequency  crosses  the  top 
edge  of  the  low  wave  number  limit  of  the  spin  wave  band. 
More  specifically,  the  FMR  point  starts  above  this  limit  at 
low  field  and  moves  into  the  band  only  at  high  field.  For 
Y1G  spheres,  as  an  example,  this  band  edge  crossover  point 
occurs  in  the  3.5  GHz  range,  depending  on  the  actual 
saturation  induction  of  the  material 

Even  today,  many  FMR  experts  do  not  know  about  this 
basic  effect,  even  though  it  provides  the  historical 
experimental  key  to  two  magnon  scattering  relaxation.  One 
possible  reason  is  that  the  predominant  FMR  work  in  recent 
years  has  been  for  thin  films.  This  band  edge  crossover 
effect  occurs  only  for  spheres  and  not  for  any  of  the  other 
basic  geometries  of  interest. 

The  experimental  manifestation  of  this  effect  is  an 
abrupt  increase  in  the  linewidth  as  the  FMR  frequency 
moves  through  the  crossover  frequency,  taken  here  as  fx . 
This  is  then  followed  by  a  slow  decrease  in  the  linewidth  at 
higher  frequencies.  This  is  a  purely  geometric  effect.  It 
occurs  for  spheres,  but  not  for  other  sample  shapes.  C.  R. 


Buffler  was  the  first  to  observe  this  effect  experimentally. 
This  direct  early  experimental  evidence  for  TMS  processes 
provided  the  key  data  that  pointed  the  way  to  further  theory 
and  experiments  on  anisotropy,  porosity,  and  surface  pit 
TMS  processes.11  As  noted  in  the  introduction,  new  results 
continue  to  emerge.  Mo  and  co-workers  have  shown,6  for 
example,  that  grain  boundary  TMS  processes  related  to  very 
short  wavelength  spin  waves  must  be  added  to  the  above 
list  of  interactions  for  a  complete  picture  of  this  important 
mechanism. 

Figure  1  illustrates  the  band  edge  crossover  effect 
described  above.  The  (a),  (b),  and  (c)  graphs  show  spin 
wave  band  diagrams  of  normalized  spin  wave  frequency  vs. 
wave  number  k  for  three  different  fields  for  sphere  shaped 
Y1G  samples,  as  indicated.  Graph  (b)  corresponds  to  the 
FMR  band  edge  crossover.  The  specific  field/ frequency 
operating  points  were  chosen  to  demonstrate  this  effect. 
The  dashed  horizontal  line  shows  the  degenerate  mode  cut 
across  the  bands  at  the  FMR  frequency  for  these  three 
cases.  The  vertical  shaded  region  across  the  low  k  edge  of 
the  band  for  k  <  ka  =  lxlO4  rad/cm  accentuates  those 
modes  that  can  give  a  strong  GG  TMS  contribution  for  Y1G 
grain  sizes  in  the  10  pm  range.  The  connection  between 
ka ,  grain  size,  and  TMS  processes  will  be  elaborated 
shortly. 

Figure  1  is  intended  to  convey  two  main  points.  First, 
as  one  increases  the  external  field,  the  well  known  upward 
shift  in  spin  wave  band  and  the  FMR  frequency  is 
accompanied  by  a  crossover  of  the  FMR  frequency  from 
outside  the  low  k  part  of  the  band  at  low  field  to 
completely  inside  the  band  at  high  field.  The  crossover 
frequency  fx  is  equal  to  (2 /3)(\y\ 1 2n)4  ?rMs  ,  where  y 
denotes  the  electron  gyromagnetic  ratio.  With  4 nMs 
taken  at  the  nominal  Y1G  value  of  1750  G  and  \y\l2n 
taken  at  the  nominal  free  electron  value  of  2.8  GHz/kG,  this 
crossover  effect  occurs  at  about  3.3  GHz.  The  specific 
diagrams  in  Fig.  1  have  been  selected  to  underscore  this 
basic  effect.  It  is  a  little  known  effect,  even  though  it  was 
crucial  to  the  initial  realization  of  the  role  of  TMS  processes 
in  microwave  ferrite  relaxation. 

The  second  point  concerns  the  interplay  between  the 
actual  coupling  to  the  degenerate  modes  and  the  crossover 
effect  that  leads  to  the  TMS  signature  in  the  linewidth  vs. 
frequency  profile  discovered  by  Buffler  In  this  context,  the 
key  is  in  the  shaded  vertical  strip  of  modes  at  the  extreme 
low  k  side  for  each  of  the  spin  wave  band  diagrams  in  Fig. 
1.  This  strip  corresponds  to  the  modes  that  are  strongly 
coupled  to  the  FMR  mode  for  coarse  grain  ferrites.  Here, 
“coarse  grain”  is  taken  to  denote  grain  sizes  in  the  range  of 
tens  of  micrometers  that  are  typical  for  polycrystalline 
ferrites  made  by  standard  sintering  or  hot  pressing  methods. 
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(a)  Low  field 


(b)  Band  edge 


(c)  High  field 


FIG  1. (Color  online)  Schematic  spin  wave  band  diagrams  of  spin  wave  frequency  fk  ,  normalized  to 
the  ferromagnetic  resonance  frequency  (FMR)  ,  as  a  function  of  the  spin  wave  number  k  for 
sphere  shaped  samples.  Graphs  (a),  (b),  and  (c)  are  for  external  applied  static  field  ( //ext )  values  of 
0.6,  1.2  and  3.6  kOe  and  the  corresponding  FMR  frequency  values,  as  indicated.  Graph  (b) 
corresponds  to  the  band  edge  crossover  point  at  =  fx  .  The  curves  are  for  yttrium  iron  garnet 
magnetic  parameters.  The  curved  lines  show  the  band  limits  for  spin  wave  propagation  angles  (0k) 
relative  to  the  field  direction  of  zero  and  /r/2,  as  indicated.  The  shaded  regions  indicate  the  full 
band.  The  horizontal  dashed  lines  show  the  cut  across  the  band  at  the  FMR  frequency  for 
fk  I  f  fmr  =  1  •  Tta  arrow  at  ka  ~  1  x  104  rad/cm  and  narrow  vertical  shaded  regions  for  k  <  ka  denote 
the  region  of  strong  two  magnon  grain-to-grain  scattering  for  a  nominal  grain  size  of  10  |im  . 


Generally,  two  magnon  scattering  has  a  strong  coupling 
between  the  driven  FMR  mode  for  k  «  0  to  degenerate  spin 
wave  modes  with  k  <  ka  »  \  l  a  or  so,  where  a  is  the  grain 
size.  For  a-  10  pm,  one  has  ka  *  104  rad/cm.  One  can 
see  that  as  the  field  is  increased  and  the  FMR  frequency 
moves  from  above  the  band  edge  point  to  inside  the  band, 
the  low  k  modes  available  for  strong  scattering  come  into 
play  quite  abruptly.  The  fact  that  the  density  of  states  for 
the  degenerate  spin  wave  modes  is  also  quite  large  at  the 
top  edge  of  the  spin  wave  band  leads  to  a  pronounced  peak 
in  the  linewidth  as  one  moves  through  the  crossover 
frequency.  The  quantitative  connection  between  the 
linewidth  vs.  frequency  data  and  the  TMS  analysis  for  GG 
scattering  will  be  covered  in  Section  IV.  Additionally,  the 
results  will  also  reveal  a  new  effect,  namely,  a  nonzero  GB 
scattering  to  high  k  spin  wave  modes.  This  effect  shows 
up  most  clearly  for  frequencies  in  the  band  regime  of  graph 
(a)  in  Fig.  1.  One  can  see  that  for  GG  scattering  alone,  the 
TMS  linewidth  in  this  regime  is  zero  because  there  are  no 
degenerate  low  k  modes. 

III.  SAMPLES  AND  MEASUREMENTS 

The  YIG  samples  were  made  by  standard  hot  isostatic 
pressing  (Hipping)  ceramic  methods.  Nazarov  et  al.  have 
provided  the  fabrication  details.3  These  authors  also 
provided  X-band  FMR  linewidth  data  that  establish  the 
density  of  the  materials  to  be  near  theoretical.  The  residual 


porosity  was  less  than  1%  and  the  half  power  10  GHz  FMR 
linewidth  was  13  Oe.  The  average  grain  size  was  8  //m . 
Spheres  were  fabricated  from  the  interior  regions  of  the 
hipped  blocks,  to  avoid  possible  problems  with  oxygen 
deficient  surface  regions.  A  similar  procedure  was  also 
applied  to  Ca-Ge  substituted  YIG.  The  specific  results  in 
this  chapter,  shown  for  nominal  2  mm  diameter  spheres, 
confirm  the  nearly  complete  elimination  of  porosity  for  the 
hipped  materials. 

The  FMR  profiles  were  measured  at  frequencies 
between  1.95  and  6  GHz  with  a  broadband  strip 
transmission  line  FMR  spectrometer  described  in  Ref.  [12]. 
Figure  2  shows  representative  data.  The  graph  shows  a 
sequence  of  profiles  of  the  FMR  absorption  derivative  vs. 
static  external  field  HexX  for  selected  frequencies,  as 
indicated.  One  sees  the  expected  shift  in  the  FMR  fields 
with  increasing  frequency.  Note  that  the  profile  is 
somewhat  narrow  at  1.9  GHz,  becomes  very  narrow  at  3 
GHz.  broadens  substantially  at  4.0  GHz,  and  then  appears 
to  narrow  somewhat  at  5.3  GHz.  The  behavior  from  3  to 
5.3  GHz  is  a  direct  consequence  of  the  band  edge  crossover 
GG  TMS  effect  from  Sec.  II.  The  substantial  jump  in 
linewidth  from  3  to  4  GHz  is  a  direct  result  of  this 
crossover.  The  fall-off  from  4  to  5.3  GHz  is  due  to  a 
gradual  drop  in  the  available  degenerate  spin  wave  states 
for  scattering.  The  relatively  broad  line  at  1.9  GHz  is  due 
to  the  unsaturated  state  of  the  sample  at  these  fields 
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Applied  field  (Oe) 

FIG.  2  Normalized  ferromagnetic  resonance  (FMR)  absorption 
derivative  vs.  field  profiles  for  the  hipped  YIG  sphere  at  selected 
FMR  pump  frequencies,  as  indicated. 

Data  of  the  sort  shown  in  Fig.  2  were  also  used  to 
establish  the  gyromagnetic  ratio  y .  The  FMR  frequency 
/fmr  f°r  sphere  shaped  samples  is  equal  to  (\y\l2n)H  . 
The  data  on  FMR  frequency  vs.  field  yielded  a  \y\I2tt 
value  of  2.765  MHz/Oe.  This  paper  will  use  Gaussian  units 
for  the  working  equations  for  the  data  analysis,  as  needed, 
and  for  the  theoretical  development.  Note  also  that  the 
44  /  ”  notation  will  denote  frequency  in  Hz  or  related  units 
In  the  theoretical  discussions  to  follow,  it  will  also  be 
appropriate  to  bring  in  the  usual  44 co”  notation  for  angular 
frequency  in  rad's. 

Linewidth  values  were  obtained  from  FMR  derivative 
profiles  similar  to  those  shown  in  Fig.  2,  but  for  a  fine  grid 
of  frequency  points  from  2  to  6  GHz.  For  each  profile,  the 
half  power  field  swept  linewidth  AH  was  taken  as  the 
difference  in  field  values  at  the  extrema  of  the  derivative 
profile  multiplied  by  V3  .  This  connection  is  strictly 
applicable  to  Lorentzian  absorption  profiles.  Integrated 
absorption  profiles  were,  in  fact,  near  Lorentzian. 

Figure  3  shows  representative  AH  vs.  results. 

These  particular  data,  as  in  Fig.  2,  are  for  the  hipped  YIG 
sphere.  The  main  portion  of  the  data  is  shown  by  solid 
symbols.  These  data  will  be  the  focus  of  the  two  magnon 
analysis  to  follow.  The  portion  of  the  data  for  <  3 
GHz  or  so  are  shown  by  blue  open  circles.  These  data  show 
a  gradual  upturn  below  3  GFiz  and  a  more  rapid  increase 
below  about  2.2  GHz.  This  response  is  related  to 
demagnetization  and  low  field  loss  effects  and  is  not  of 
direct  interest  here.  The  vertical  dashed  line  at  *3.6 
GHz  marks  the  band  edge  crossover  frequency.  This  point 
is  shifted  from  the  3.3  GHz  value  shown  in  Fig.  1  because 
of  the  larger  nominal  4 nMs  value  of  1 950  G  for  the  hipped 
YIG  material,  as  obtained  from  fits  to  be  considered  shortly. 
The  schematic  low  field  and  high  field  region  spin  wave 


Band  edge  crossover  frequency 


FMR  frequency  f  FMR  (GHz) 


FIG.  3  (Color  online)  Ferromagnetic  resonance  (FMR)  linewidth 
AH  as  a  function  of  the  FMR  frequency  /FMR  for  the  hipped  YIG 
sphere  The  open  blue  symbols  show  data  for  the  low  field  regime 
for  which  the  sample  is  not  saturated.  The  vertical  dashed  line 
marks  the  calculated  band  edge  crossover  frequency  fx  .  The 
insets  show  schematic  spin  wave  frequency  vs.  wave  number 
diagrams  for  the  low  field  below  crossover  ( /  <  fx  )  and  high 
field  above  crossover  f  >  fx  regimes. 

band  diagrams  for  the  regions  below  and  above  the 
crossover  point  track  the  (a)  and  (c)  diagrams  from  Fig.  1 . 

The  solid  circle  data  in  Fig.  3  demonstrate  the  two 
magnon  band  edge  crossover  effect  discussed  in  Sec.  II. 
There  is  a  rapid  increase  in  linewidth  as  one  approaches  the 
crossover  frequency  from  below,  followed  by  a  distinct 
peak  and  a  gradual  leveling  off  at  high  frequency.  It  is 
notable  that  there  is  sizeable  linewidth  increase  even  for 
frequencies  below  the  crossover  frequency.  The  coarse 
grain  scattering  scenario  from  Sec.  II  and  the  exclusive  two 
magnon  coupling  to  extremely  low  k  modes  alone  would 
give  a  sharp  cutoff  at  the  crossover  frequency  fx  rather 
than  the  response  shown  by  the  data.  As  will  be  discussed 
shortly,  the  two  magnon  linewidth  “leakage”  into  the 
frequency  regime  below  fx  is  a  direct  result  of  grain 
boundary  TMS  processes.  It  also  turns  out  that  GB 
scattering  also  plays  a  significant  role  in  the  frequency 
regime  above  crossover.  This  is  alluded  to  in  the  last  part 
of  the  discussion  of  Fig.  1  at  the  end  of  Sec.  II.  These 
points  will  be  made  quantitative  in  Secs.  IV  and  V. 

IV  TWO  MAGNON  SCATTERING 
RELAXATION  AND  LINEWIDTH  CONNECTIONS 

As  noted  in  Section  II,  there  is  extensive  literature  on 
the  two  magnon  scattering  relaxation  theory  in  magnetic 
systems  in  general.8’11  15  Reference  4  casts  the  general 
two  magnon  scattering  relaxation  rate  in  the  form 
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^TMS  =  j  <?(<y-^)C(k)A(M<*k-  O 

k-space 

The  parameter  coM  =  j ;  1 4 nMs  expresses  the  saturation 
induction  in  frequency  units.  The  delta  function  S(co-cok) 
constrains  the  scattering  to  degenerate  spin  wave  modes 
with  a  frequency  cok  equal  to  the  pump  frequency  co  .  The 
C(k)  factor  defines  the  coupling  strength  between  the 
uniform  mode  and  particular  degenerate  spin  wave  modes 
at  some  wave  vector  k  .  The  A(k,tf>)  factor  accounts  for 
the  ellipticity  of  the  spin  wave  mode.  For  sphere  shaped 
samples,  the  connection  between  /7tms  A//  amounts 
to  a  simple  conversion,  according  to  A II  =  2rjTm  f\y\. 

As  discussed  in  earlier  sections,  two  distinct  TMS 
processes  turn  out  to  be  important  for  polycrystalline 
ferrites,  one  related  to  GG  scattering  and  to  GB  scattering. 
These  processes  can  be  taken  to  occur  in  parallel  In  order 
to  deal  with  the  differences  in  the  dependences  of  the 
coupling  terms  for  the  two  processes,  it  proves  useful  to 
write  Eq.  (1)  in  a  slightly  different  form.  Keep  in  mind  that 
for  GG  scattering,  it  is  the  grain  size  a  that  controls  the 
scattering  and  the  coupling  is  to  low  wave  number  modes 
with  k  <  ka  alone.  Recall  the  discussion  related  to  Fig.  1 
and  coarse  grain  scattering.  For  GB  scattering,  it  is  the 
much  narrower  sub-micron  size  grain  boundary  that 
controls  the  scattering.  This  process  is  definitely  not  coarse 
grain  scattering  but  extends  out  to  much  larger  k  -  values. 

A  useful  starting  point  for  the  separation  of  the  different 
k  -  dependences  for  GG  and  GB  scattering  is  to  recast  the 
formal  k-space  integral  in  Eq.  (1)  in  terms  of  the  specific 
integrals  over  wave  number  k  and  the  polar  and  azimuthal 
spin  wave  propagation  angles  0^  and  ^  ,  respectively. 


The  new  g(k)  function  is  used  here  to  separate  out  the 
k  -  dependent  part  of  the  degenerate  spin  wave  coupling  for 
GG  scattering.  This  factor  is  taken  in  the  form  first  used  by 
Schlomann'  and  recently  re-derived  by  hamiltonian 
methods  by  Krivosik  et  al.  in  Ref.  14 


«(*)= 


4  kW 

»(i  +  *V)2’ 


(4) 


In  the  above,  H A  denotes  an  effective  cubic  anisotropy 
field  2K\IMS,  where  Ky  is  the  first  order  cubic 
magnetocrystal  line  anisotropy  energy  constant.  In  the 
qualitative  discussion  above,  a  was  taken  as  the  grain  size. 
In  a  more  general  sense,  a  may  be  taken  as  a  mean 
inhomogeneity  size.  One  can  see  that  the  g(£)  function 
peaks  at  k  =  ka  =  \  l a  .  Recall,  however,  that  the  k2  factor 
in  the  numerator  derives  from  the  same  factor  in  the 
k  -  space  differential  volume  4 nk2dk  .  Without  this  factor, 
the  g(k)  corresponds  to  the  Fourier  transform  of  a 
fluctuating  field  with  a  correlation  length  equal  to  the  mean 
grain  size. 

The  spin  wave  ellipticity  factor  for  GG  scattering, 
labeled  here  as  AGG(k,^)  was  also  developed  first  by 
Schlomann  in  Ref.  9  and  is  written  here  as 

,  x  A(k,H) 

Ago  (k,fl>)  =  1  +  19 — — — (5) 

CO 


where 


A (k,  H )  =  \y\  H,  +  \y\ Dk2  +  2ttMs  sin2  6k .  (6) 


2  00  1 

7tms  =^j  jdk  Jdcos0*[F(k)A(k,<»)  S(o)-o)k)]. 
87r  0  0 

(2) 

The  function  F(k)  =  4nk2C(k)  folds  in  the  numerical  4n 
factor  from  the  angular  integrals  and  the  explicit  k 2  factor 
from  the  A -integral  into  the  new  form  of  the  coupling.  As 
the  discussion  below  will  demonstrate,  the  F^(k)  form 
makes  it  possible  to  gain  an  intuitive  understanding  of  the 
GG  and  GB  contributions  to  the  scattering. 

First  consider  the  working  equations  for  GG  scattering. 
In  this  case,  one  can  write  ^(k)  in  the  form 


^GG  00  -  ^GG  (^)  " 


8  7T4 

105 


H. 


4ttM 


g(k)- 


S  / 


(3) 


The  internal  field  Hi  in  Eq.  (6)  is  equal  to  II  -(4 ttMs  /3) 
for  sphere  shaped  samples.  The  D  parameter  takes 
exchange  into  account.  The  nominal  D  -value  for  YIG 
materials  is  5.2xl0~9  Oe  cm2.4  The  numerical  factors  in 
Eqs  (3)  and  (5)  follow  the  Schlomann  form.  They  are  a 
result  of  a  systematic  averaging  over  the  randomly  oriented 
cubic  crystallites  for  the  YIG  system.  Further  details  on  the 
averaging  mechanics  are  given  by  McMichael  and  Krivosik 
in  Ref.  13. 

For  GB  scattering,  the  F(k)  and  A(k,rj)  functions 
take  a  somewhat  different  form.  The  full  theory  is 
summarized  in  Refs.  4  and  6,  and  developed  in  detail  in 
Ref.  15.  The  coupling  and  polarization  terms  are  different 
from  those  for  the  grain-grain  anisotropy  based  scattering. 
The  formulae  given  below  follow  the  nomenclature  of  Ref. 
4.  The  F(k)  for  GB  scattering  now  takes  the  form 
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^gb  (k)=Fw(k)  =  l2x*ll„k1g(k).  (7) 

Following  Ref.  4,  lgt  denotes  an  effective  grain  boundary 
thickness  parameter.  From  physical  considerations,  lgb 
scales  with  the  ratio  of  surface  anisotropy  energy  to  the 
magnetostatic  self  energy  and  is  on  order  of  tens  of 
nanometers  for  polycrystalline  ferrites.  Fits  to  the  data  in 
the  next  section  will  yield  nominal  values  of  9  nm  for  the 
pure  YIG  and  16  nm  for  the  Ca-Ge  substituted  Y1G 
samples.  Note  the  additional  A 2  factor  in  the  GB  coupling. 
With  this  modification,  the  A;  -  dependence  of  the  scattering 
defined  through  A2/(A)  sees  major  modifications.  Now 
the  scattering  is  essentially  zero  for  k  «  ka  ,  begins  to  grow 
for  A  a  ka ,  and  levels  off  to  some  maximum  value  for 
k  »  ka .  This  is  precisely  the  response  expected  for  GB 
scattering.  One  can  now  have  strong  coupling  to  very  short 
wavelength  spin  waves  in  the  extremely  high  k  part  of  the 
spin  wave  manifold. 

Also  following  Ref.  4,  the  A(k,6>)  spin  wave 
polarization  term  for  GB  scattering  is  now  written  as 

/  v  A  k,//) 

Agb  (k,&>)  =  1  +  V  (8) 

0) 

The  different  numerical  factors  in  Eqs.  (6)  and  (8)  arise 
from  the  different  symmetries  in  the  GG  and  GB  scattering 
problems.  As  noted  above,  GG  scattering  involves  an 
averaging  over  the  grains  in  the  3D  system  For  GB 
scattering,  in  contrast,  one  is  dealing  with  an  interface  effect 
that  is  treated  in  a  simplified  2D  model.1' 

Figure  4  shows  the  same  bulk  spin  wave  band  diagrams 
of  frequency  f  k  vs.  wave  number  k  as  given  in  graphs  (a) 
and  (c)  of  Fig.  1,  except  that  the  A -axis  is  given  in 
logarithmic  scale.  Also  as  in  Fig.  1,  horizontal  dashed  red 
lines  are  used  to  show  k-  range  of  the  degenerate  mode  at 
the  indicated  FMR  frequency  in  each  case.  The  solid  blue 
and  dashed  green  curves  in  both  graphs  show  plots  of  the 
F(k)  functions  developed  above  for  GG  and  GB 
scattering,  respectively.  The  curves  are  normalized  to 
maximum  values  and  denoted  as  Fnorm(A) .  Recall  that  the 
peak  GG  coupling  occurs  for  k  =  ka  =  \/a.  The  GB 
coupling  shows  no  peak,  but  simply  levels  off  at  its 
maximum  value  at  very  high  A  -  values. 

The  overall  point  of  Fig.  4  is  to  demonstrate  the  very 
different  coupling  conditions  for  fields  below  the  band 
edge,  as  in  (a),  and  fields  above  the  band  edge,  as  in  (b). 
For  the  below-the-band-edge  situation  in  (a),  one  can  see 
that  the  degenerate  mode  line  cuts  the  spin  wave  band  in  the 
range  of  very  high  A  values  for  which  the  GB  coupling  is 
strong  and  that  there  are  no  degenerate  low  k  modes.  This 
means  that  the  two  magnon  processes  are  limited  to  GB 


10°  io'  102  103  10*  105 


Spin  wave  number  k  (rad/cm) 

FIG  4.  (Color  online)  Each  graph  shows  (i)  spin  wave  frequency 
fk  vs.  wave  number  k  for  a  sphere  shaped  yttrium  iron  garnet 
sample  with  materials  parameters  that  match  the  experiment  and 
(ii)  the  normalized  two  magnon  scattering  coupling  coefficient 
F"0""  vs.  k  for  grain-to-grain  (GG)  and  grain  boundary  (GB)  two 
magnon  scattering  (TMS),  as  indicated.  The  curves  in  (a)  and  (b) 
were  obtained  with  external  field  //ext  -values  of  0.6  and  3.6  kOe, 
respectively.  The  spin  wave  band  nomenclature  and  labels  are  the 
same  as  for  Fig.  1.  The  solid  vertical  line  at 
A  =  ka  =  1  la  =  1 .25  x  1 03  rad/cm  ,  where  a  is  the  nominal  grain 
size,  marks  the  transition  k  -  value  for  both  TMS  processes. 

scattering.  As  will  be  shown  shortly,  GB  scattering  turns 
out  to  provide  a  reasonable  quantitative  model  for  the 
observed  linewidth  tail  for  frequencies  below  the  band  edge 
crossover  frequency  f x  . 

Contrast  this  with  the  situation  in  Fig  4  (b).  Here  one 
can  see  that  the  degenerate  mode  line  cuts  the  spin  wave 
band  over  the  full  range  of  wave  number  A  values  from  the 
low  limit  at  A  «  0  out  to  the  very  high  values  at  the  edge  of 
the  spin  wave  band  for  A»  105  -  106  rad/cm.  These 
degenerate  modes  now  span  regions  of  A -space  for  which 
the  F(A)  functions  for  both  GG  and  GB  scattering  indicate 
strong  coupling.  The  quantitative  fits  of  the  theory  to  the 
data  given  below  will  show  that  both  processes  contribute 
substantially  to  the  two  magnon  linewidth  for  frequencies 
above  the  band  edge  crossover  frequency  fx  . 
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V  LINEWIDTH  ANALYSIS 

Graphs  (a)  and  (b)  in  Fig.  5  show,  respectively,  the  full 
ensemble  of  linewidth  vs.  frequency  data  for  both  the  pure 
YIG  and  substituted  YIG  samples.  The  open  circles  in  (a) 
show  the  same  data  as  in  Fig.  3.  The  solid  circles  show 
extended  frequency  data  from  Ref.  3.  The  error  bars  for  all 
data  points  are  in  the  range  of  ±0.5  Oe.  The  panels  also 
show  theoretical  fits  based  on  a  simple  model  of  GG  and 
GB  TMS  processes  that  operate  in  parallel.  The  fits  do  not 
include  either  a  Landau-Lifshitz  linewidth  term  or 
inhomogeneity  broadening  effects.  The  various  curves 
show  computed  results  for  GG  scattering,  GB  scattering, 
and  the  sum  of  these  two  processes,  as  indicated.  The 
theoretical  curves  were  obtained  for  a  common  grain  size  a 
of  8  /im  ,  and  with  4 nM s ,  // A  ,  l  b  ,  and  \y\!2n  values  of 
1952  G,  40  Oe,  9  nm,  and  2.765  MHz/Oe,  respectively,  for 
the  pure  YIG  and  990  G,  20  Oe,  16.2  nm,  and  2.76 
MHz/Oe,  respectively,  for  the  substituted  YIG.  These 
materials  parameters  are  consistent  with  the  ferrite 
literature.1  1S  The  critical  fit  parameters  were  H A  and  lgb  . 
Following  Ref.  4,  an  effective  grain  boundary  thickness 
parameter  lgb  in  the  nm  range  is  also  reasonable  for  the 
hipped  ferrite  microstructure. 

Both  graphs  show  very  prominent  Buffler  peaks.  In  line 
with  the  extended  discussion  in  Sec.  Ill,  these  peaks  provide 
the  basic  signature  of  GG  scattering.  The  prominence  of 
the  Buffler  peaks  is  due  to  the  use  of  hipped  samples  with 
near  theoretical  density.  Porosity  effects  in  the  original 
Buffler  samples  resulted  in  significantly  larger  linewidths 


and  a  broader  peak.  Overall,  the  linewidths  in  (b)  for 
sample  B  are  even  lower  than  those  in  (a)  because  of  the 
use  of  Ca-Ge  substituted  YIG  materials  with  a  lower 
anisotropy.  The  GG  based  two  magnon  scattering  scales 
with  the  anisotropy.16,17 

The  physical  origins  of  the  linewidth  vs.  frequency 
response  manifested  by  the  low  frequency  data  in  (a)  was 
discussed  in  connection  with  Fig.  3.  From  Fig.  5,  one  sees 
that  both  samples  show  essentially  the  same  response.  Each 
data  set  shows  (1)  a  rapid  increase  in  linewidth  at  very  low 
frequencies  due  to  sample  demagnetization,  (2)  a  local 
minimum  with  a  nonzero  linewidth  at  frequencies  just 
below  the  Buffler  peak,  (3)  a  rapid  increase  that  forms  the 
Buffler  peak  as  one  moves  through  the  band  edge 
crossover,  and  (4)  a  gradual  decrease  in  linewidth  as  one 
moves  to  frequencies  well  above  the  Buffler  peak.  These 
effects  occur  for  somewhat  lower  frequencies  in  (b)  due  to 
the  reduced  magnetization  for  the  Ca-Ge-YIG. 

The  new  points  to  be  considered  here,  relative  to  Fig.  3, 
concern  the  actual  fits  to  the  data  from  the  combined  GG 
and  GB  TMS  models.  These  fits  tell  the  quantitative  story 
of  the  TMS  processes  outlined  in  general  terms  in  Secs.  Ill 
and  IV.  The  fits  demonstrate  three  points.  First,  it  is  clear 
that  one  can  realize  a  remarkably  good  fit  to  the  data  for 
two  different  materials  based  on  TMS  processes.  The  fact 
that  the  theoretical  curves  show  a  divergence  at  the  band 
edge  crossover  frequency  is  of  no  consequence  here.  This 
divergence  is  due  to  the  neglect  of  the  nonzero  relaxation 
rate  for  the  product  spin  waves  from  the  TMS  processes. 
SchlomanrT  has  shown  that  the  inclusion  of  a  nonzero  spin 


FIG.  5  (Color  online)  Graphs  (a)  and  (b)  show  the  ferromagnetic  resonance  (FMR)  linewidth  AH  as  a  function  of  the  FMR  frequency  /FMR 
for  sphere  shaped  samples  of  hipped  yttrium  iron  garnet  (YIG)  and  Ca-Ge  substituted  YIG,  respectively.  The  open  circles  show  data  from  the 
broadband  strip  transmission  line  measurements.  The  solid  circles  in  (a)  show  data  on  the  same  sample  from  Nazarov  et  al ?  The  dashed  and 
dotted  curves  show  fitted  AH  vs.  /FMR  results  from  model  calculations  for  grain  boundary  (GB)  and  grain-to-grain  (GG)  two  magnon 
scattering  processes,  respectively.  The  solid  curves  show  the  sum  of  the  GB  and  GG  results. 
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waves  linewidth  in  TMS  calculations  eliminate  such 
divergences.  One  obtains,  thereby,  theoretical  curves  that 
match  effective  linewidth  data  quite  nicely,  for  example.1’ 

Second,  the  GB  TMS  linewidth  fits  for  the  frequency 
region  just  below  the  Buffler  peak  band  confirm  the 
expectation  from  Sec.  IV  that  grain  boundary  scattering 
processes  must  play  a  role  in  this  region.  In  fact,  the 
theoretical  response  also  matches  this  portion  of  the  data 
quite  nicely. 

The  third  point  comprises  a  somewhat  unexpected 
result.  The  realization  of  data  over  a  wide  range  of 
frequencies  on  both  sides  of  the  Buffler  peak,  in 
combination  with  calculations  based  on  a  two  component 
GG/GB  TMS  process;  show  that  both  processes  make 
significant  contributions  to  the  linewidth.  The  fits  in  Ref.  3 
were  done  for  GG  scattering  alone  and  appeared  to  match 
the  data.  The  problem  in  this  reference,  however,  was  in 
the  availability  of  limited  high  frequency  data  only  and  no 
experimental  points  close  to  or  below  the  Buffler  peak.  In 
the  present  work,  it  is  these  low  frequency  data  points  that 
provide  the  clue  to  the  role  of  GB  scattering.  Once  the 
required  GB  scattering  parameters  are  included  to  fit  the 
data  for  f  <  fx  ,  a  necessary  GB  contribution  propagates 
into  the  f  >  fx  frequency  regime.  The  fits  show  that  in 
the  /  >  fx  region,  one  has  roughly  equal  GG  and  GB  two 
magnon  scattering  contributions  to  the  linewidth. 

With  two  magnon  processes  properly  taken  into 
account,  it  is  entirely  reasonable  to  ignore  any  possible 
intrinsic  linewidth  terms  in  doing  the  fits  shown  in  Fig  5. 
Effective  linewidth  measurements  have  shown  that  such 
contributions  are  below  one  Oersted  over  the  frequency 
range  shown  here.19 

VI  SUMMARY  AND  CONCLUSION 

The  early  linewidth  results  for  polycrystalline  ferrites  in 
the  region  of  the  Buffler  peak  provided  the  initial  evidence 
for  the  role  of  two  magnon  scattering  relaxation  in 
ferromagnetic  resonance.  These  classic  results  have  been 
now  extended  and  made  quantitative,  and  several  new 
aspects  of  TMS  relaxation  in  polycrystalline  ferrites  have 
been  revealed.  The  new  data  on  ultradense  pure  and 
substituted  YIG  materials  with  nominal  10  micron  grain 
sizes  show  an  especially  sharp  Buffler  peak  that  derives 
from  the  band  edge  crossover  effect  initially  reported  in 
Ref.  7.  Due  in  part  to  recent  new  developments  in  two 
magnon  relaxation  models  for  both  grain-to-grain 
anisotropy  scattering  (GG  scattering)  and  grain  boundary 
(GB)  scattering,  it  has  been  possible  to  obtain  quantitative 
TMS  fits  to  the  data.  These  fits  demonstrate  the  expected 
role  of  GG  scattering  in  the  relaxation,  just  as  suggested  by 
Buffler.  The  fits  also  reveal  the  somewhat  unexpected  role 
of  GB  scattering  as  a  contributor  to  the  TMS  linewidth  in 


polycrystalline  ferrites.  (1)  It  is  found  that  GB  scattering 
plays  a  role  in  the  frequency  regime  below  the  Buffler  peak 
where  simple  two  magnon  low-  k  scattering  models  give  no 
TMS  contribution.  The  nominal  GB  TMS  linewidth  in  this 
regime  is  about  5  Oe.  (2)  It  is  also  found  that  GB  scattering 
contributes  about  half  of  the  overall  TMS  linewidth  in  the 
frequency  regime  above  the  peak  where  GG  scattering  was 
previously  though  to  be  dominant. 

ACKNOWLEDGEMENTS 

This  work  was  sponsored  in  part  by  the  Office  of  Naval 
Research  Grant  Nos.  N000 14-07-1 -0597  and  N000 14-08-1- 
1050  (subgrant  from  Virginia  Commonwealth  University), 
U.S.  Army  Research  Office  Grant  No.  W91  INF-04- 1-0247 
(MUR1),  the  Information  Storage  Industry  Consortium 
extremely  high  density  recording  program,  and  Seagate 
Technology.  The  hipped  pure  YIG  sample  was  kindly 
provided  through  a  collaboration  with  G.  M.  Argentina 
(Pacific  Ceramics,  Inc.)  and  J.  J.  Green  and  H.  J.  Van  Hook 
(private  consultants,  Lexington,  Massachusetts).  The 
hipped  substituted  YIG  sample  was  kindly  provided  by  TCI 
Ceramics,  Bethlehem,  Pennsylvania. 

Footnote: 

*:  2  Now  at  CIS,  Tyco  Electronics,  Reno  Nevada  89521, 
U.S. A. 


REFERENCES 

]H  J  Van  Hook  and  C.  B.  Willingham,  Adv.  Ceram  15,  1637 
(1984). 

2H.  V.  Atkinson  and  S.  D.  A,  Metall.  Mater.  Trans.  A  31  A,  2981 
(2000). 

3A.  V.  Nazarov,  D.  Menard,  J.  J.  Green,  C.  E.  Patton,  G.  M. 

Argentina,  and  H.  J.  Van  Hook,  J.  Appl.  Phys.  94,  7227  (2003). 
4N.  Mo,  J.  J.  Green,  P.  Krivosik,  and  C.  E.  Patton,  J.  Appl  Phys 
101,023914  (2007). 

5S.  S.  Kalarickal,  P  Krivosik,  J  Das,  K.  S.  Kim.  and  C.  E  Patton, 
Phys.  Rev.  B  77,  054427  (2008). 

6N.  Mo,  J.  Hohlfeld.  M.  u  Islam,  C.  S.  Brown,  E  Girt,  P. 
Krivosik,  W.  Tong,  A.  Rebei,  and  C.  E.  Patton,  Appl.  Phys.  Lett. 
92,  022506  (2008). 

7C.  R.  Buffler,  J.  Appl  Phys  30  I72S  (1959). 

8A.  M.  Clogston,  H  Suhl,  L.  R.  Walker,  and  P  W  Anderson,  J 
Phys.  Chem.  Solids  1,  129  (1956). 

9E.  Schloemann,  J.  Phys.  Chem.  Solids  6,  242  (1958). 

,0E.  Schloemann,  Phys.  Rev.  182,  632  (1969). 
nM.  Sparks,  Ferromagnetic  Relaxation  Theory  (McGraw-Hill, 
New  York,  1964). 

,2S.  S.  Kalarickal,  P.  Krivosik,  M.  Wu,  C.  E.  Patton,  M.  L. 
Schneider,  P.  Kabos,  T.  J.  Silva,  and  J.  P.  Nibarger,  J.  Appl  Phys 
99,  093909  (2006). 


Kalarickal  et  al. ,  Phys.  Rev.  B  submit.  FMR  in  hipped  YIG  spheres  TMS  revised,  BZ10690 


Page  9 


13R  D.  McMichael  and  P.  Krivosik,  IEEE  Trans.  Magn.  40,  2 
(2004) 

I4P.  Krivosik,  N.  Mo,  S.  S.  Kalarickal,  and  C.  E.  Patton,  J.  Appl 
Phys  101,  083901  (2007). 

I5N.  Mo,  Ph.D.  thesis,  Colorado  State  University,  2006. 


16H  J  Van  Hook,  J.  J.  Green,  F.  Euler,  and  E  R  Czerlinsky,  J. 
Appl  Phys  39,  730  (1968). 

17C.  E.  Patton  and  H  J.  Van  Hook,  J  Appl.  Phys  43,  2872  (1972). 
l8C.  E.  Patton,  Phys.  Rev.  179,  352  (1969). 

19N.  Mo,  Y.  Y.  Song,  C.  E.  Patton,  J.  Appl.  Phys.  97,  093901 
(2005) 


Ordonez-Romero,  et  al.,  BLS  observation  of  splitting  and  confluence  spin  waves,  ver  66,  02  Feb  2009  Page  1 
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garnet  films  -  Wave  vector  selective  Brillouin  light  scattering  measurements  and  analysis. 

Cesar  L  Ordonez-Romero,1,2  Boris  A,  Kalinikos,1’3  Pavol  Krivosik,1  4  Wei  Tong,1 
Pavel  Kabos,5  and  Carl  E.  Patton1 

1  Department  of  Physics ,  Colorado  State  University ,  Fort  Collins ,  Colorado ,  USA 

2  Universidad  Nacional  Autonoma  de  Mexico,  Ciudad  Univer  sit  aria,  Mexico  D.F,  Mexico 

3  St.  Petersburg  Electrotechnical  University ,  197376  St  Petersburg,  Russia 

4  Slovak  University  of  Technology,  81219  Bratislava,  Slovak  Republic. 

5  National  Institute  of  Standards  and  Technology,  Boulder ,  Colorado,  USA 

(Distribution  version  of  manuscript  for  submission  for  publication  in  Phys.  Rev.  B,  target  submission  date  9  Feb  2009) 

Brillouin  light  scattering  (BLS)  has  been  used  to  observe  and  confirm  the  existence  of  nonlinear  three  magnon 
splitting  and  confluence  processes  for  propagating  spin  waves  in  the  magnetostatic  backward  volume  wave 
(MSBVW)  configuration.  Wave-vector  and  frequency  selective  BLS  techniques  were  also  used  to  provide  a 
quantitative  map  of  the  wave  vector  make-up  for  the  parametrically  excited  half-frequency  dipole-exchange 
spin  wave  (DESW)  split  magnons  and  the  confluence  magnons  that  result  from  the  recombination  of  these 
DESW  modes.  The  experimental  wave  vector  maps  for  the  product  splitting  and  confluence  magnons  matched 
nicely  with  those  expected  from  spin-wave  theory.  The  data  were  obtained  with  (I)  a  strip  line 
excitation/detection  transducer  structure,  (2)  forward  scattering  BLS  optics,  (3)  a  fixed  magnetic  field  of  352  Oe 
applied  along  the  propagation  direction,  (4)  pumping  frequencies  from  2.5  down  to  2.1  GHz,  (5)  and  cw  input 
powers  from  200  pW  to  6  mW.  The  wave-vector  selective  measurements  utilized  variable  diameter  circular 
diaphragms,  rotatable  slit  apertures,  and  circular  light  blocks  to  access  spin  waves  with  wave  numbers  from 
about  100  to  3.6  x  1 04  rad/cm  and  the  full  360°  range  of  propagation  angles. 

PACS  numbers:  05.45. -a,  75.30.Ds,  76.50.+g,  78.35.+C,  85.70.Ge 


I.  INTRODUCTION 

During  the  last  several  decades,  extremely  low  loss 
magnetic  thin  films,  and  yttrium  iron  garnet  (YIG)  films  in 
particular,  have  been  used  to  study  a  wide  variety  of  linear  and 
nonlinear  spin-wave  phenomena  at  microwave  frequencies.  In 
addition  to  the  relevance  of  these  effects  to  problems  in 
magnetodynamics,  the  results  are  also  of  fundamental 
importance  for  the  understanding  of  wave  behavior  in 
dispersive  dissipative  media  in  general.  Currently  there  are 
three  main  areas  of  interest  in  the  study  of  nonlinear 
magnetization  dynamics  in  magnetic  films,  (i)  microwave- 
driven  parametric  spin-wave  processes,  often  termed  Suhl 
instabilities,1'3  (ii)  modulational  instability  and  stationary  spin- 
wave  excitations  such  as  envelope  solitons,4'8  and  (iii)  chaotic 
nonlinear  dynamics.  All  three  of  these  classes  of  nonlinear 
magnetodynamics  effects  involve  nonlinear  spin-wave 
interactions  through  so-called  three  or  four  wave  processes. 

This  paper  focuses  on  item  (i),  and  in  particular,  three 
wave  splitting  and  confluence  processes  and  the  elucidation  of 
these  interactions  by  wave-vector  selective  Brillouin  light 
scattering  (BLS).  In  2000,  Synogach  et  al.  1  were  the  first  to 
demonstrate  the  role  of  parametric  three  wave  processes  in  the 
creation  of  ultra  short  microwave  pulses  in  YIG  film  spin-wave 
propagation  structures.  The  data,  while  convincing,  were 
comprised  of  microwave  results  only.  In  2003,  Mathieu  et  al.  2 
were  able  to  use  time  and  space  resolved  BLS  techniques  to 
obtain  direct  images  of  the  propagating  half-frequency  spin 
waves  produced  by  the  three  wave  splitting  interactions. 
Interestingly,  these  images  revealed  very  low  velocity  wedge 
shaped  wave  fronts  with  the  parametric  spin  waves 


concentrated  along  the  leading  edge.  These  results  had  one 
missing  link,  namely,  no  direct  evidence  for  the  companion 
confluence  process  by  which  the  half-frequency  magnons 
recombine  to  make  the  narrow  pulse  magnons  at  the  initial 
signal  frequency. 

This  missing  link  is  addressed  here.  Continuous  wave 
microwaves  were  used  to  excite  spin  waves  in  a  YIG  film 
magnetostatic  backward  volume  wave  (MSBVW)  propagation 
structure.  Carefully  measured  BLS  spectra  for  different  input 
power  levels  revealed  both  the  half-frequency  peak  that  is  the 
signature  of  the  parametric  magnons  from  the  three  wave 
splitting  as  well  as  a  prominent  signal  that  was  only  slightly 
downshifted  in  frequency  from  the  initial  pumping  frequency. 
Detailed  wave-vector  selective  BLS  measurements  served  to 
identify  this  signal  as  due  to  a  confluence  of  the  parametric 
magnons.  This  represents  the  first  direct  observation  of  three 
wave  confluence  processes  in  magnetic  systems.  It  also 
represents  the  first  direct  observation  of  splitting  for  MSBVW 
spin  waves. 

Section  II  describes  the  basic  BLS/microwave  experiment 
as  well  as  operational  details  on  wave  vector  selective 
measurements.  Section  III  presents  representative  BLS  data 
that  comprise  the  experimental  evidence  for  the  three  wave 
splitting  and  confluence  processes  and  makes  qualitative 
connections  with  the  governing  dipole  exchange  spin  wave 
(DESW)  frequency-wave  number  dispersion  relations.  Section 
IV  presents  additional  wave-vector  selective  BLS  data  that 
serve  to  identify  the  explicit  nature  of  the  splitting  and 
confluence  magnon  products  from  the  pumped  MSBVW  spin 
waves.  Section  V  provides  brief  comments  on  the  effect  of 
small  changes  in  the  pumping  frequency  on  the  magnon 
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splitting  and  confluence  processes  elucidated  in  the  earlier 
sections.  Section  VI  provides  a  summary  and  conclusions. 

II.  MICROWAVE  AND  BRILLOUIN  LIGHT 
SCATTERING  SET-UP 

The  key  microwave  component  of  the  set-up  was  the  spin- 
wave  propagation  structure.  This  was  comprised  of  a  long  and 
narrow  YIG  Film  strip  in  a  standard  MSBVW  spin-wave  (SW) 
delay  line  configuration.  The  strip  was  20  mm  long,  2  mm 
wide,  and  6.3  pm  thick.  The  low  loss  YIG  Film  was  grown  by 
standard  liquid  phase  epitaxy  techniques.  The  5  GHz 
ferromagnetic  resonance  half  power  line  width  was  about 
0.5  Oe.  A  pair  of  50  pm-wide  and  2  mm-long  microstrip 
transducers  with  a  separation  of  7  mm  was  used  to  excite  and 
detect  the  propagating  spin  waves.  For  propagation  in  the 
MSBVW  configuration,  the  static  external  Field  H  was  applied 
parallel  to  the  long  dimension  of  the  strip  and  the  SW 
propagation  direction.  The  Field  was  set  at  a  nominal  value  of 
352  Oe  for  all  of  the  measurements.  The  cw  microwave  input 
was  at  a  pumping  frequency  fp  of  2.5  GHz  for  most  of  the 
measurements.  For  the  frequency  effect  study  in  Sec.  V, 
additional  fp -values  from  2.4  to  2.1  GHz  were  also  used.  For 
the  BLS  measurements  shown  below,  the  nominal  input 
microwave  power  level  at  the  launch  transducer  ranged  from 
about  200  pW  to  6  mW.  Some  of  the  preliminary  microwave 
measurements  were  done  at  powers  down  to  one  microwatt  or 
so.  These  data  indicated  a  Suhl  threshold  power  level  in  the 
range  of  1  pW  or  so. 

The  above  Field-frequency  operating  point  was  chosen 
subject  to  two  criteria.  The  First  is  the  realization  of  a  strong 
SW  signal  at  the  output  transducer.  2.5  GHz  sits  at  the 
maximum  of  the  MSBVW  amplitude  vs.  frequency 
transmission  curve  for  the  structure  at  352  Oe.  This  operating 
point  corresponded  to  a  wave  number  k  for  the  pumped 
magnons  on  the  lowest  order  branch  MSBVW  dipole  exchange 
spin  wave  (DESW)  of  about  100  rad/cm.  The  second  is  the 
need  to  ensure  that  three  wave  splitting  processes  are  allowed. 
In  practical  terms,  this  means  that  spin  waves  at  fp! 2  are 
supported  by  the  propagation  structure.2  At  H  =  352  Oe,  the 
bottom  of  the  MSBVW  dipole-exchange  spin  wave  band  sits  at 
about  1 . 1  GHz,  well  below  fp  /  2 . 

The  BLS  measurements  were  performed  using  a  multipass 
tandem  Fabry-Perot  interferometer  (TFP)  in  a  forward 
scattering  configuration,  with  and  without  wave-vector 
selectivity.  The  basic  forward  scattering  approach  to  the 
analysis  of  pumped  spin-wave  excitations  in  YIG  films  was 
first  developed  by  Srinivasan  and  co-workers."  12  Wave- 
vector  selective  BLS  analysis  techniques  for  spin  waves  were 
first  developed  by  Wilber  and  co-workers.13'15 

The  overall  measurement  approach  used  here  follows 
Mathieu  et  al2  For  the  current  experiments,  a  514.5  nm 
wavelength  linearly  polarized  argon  ion  laser  light  was  focused 
onto  the  YIG  film  by  a  12  cm  focal  length  lens  at  normal 
incidence  with  10  mW  incident  power.  The  directly 
transmitted  and  scattered  light  was  collected  by  a  50  mm  focal 
length,  and  30  mm  diameter  FI. 4  camera  lens.  A  polarizing 
filter  set  perpendicular  to  the  incident  light  polarization  was 
used  to  reduce  substantially  the  intensity  of  the  directly 


transmitted  beam  Magnon  scattered  light  generally 
experiences  a  90°  rotation  in  polarization  from  the  incident 
beam. 

The  wave-vector  selective  BLS  measurements  were  done 
in  two  ways.  First,  different  sized  on-axis  circular  apertures 
formed  from  a  variable  diaphragm  behind  the  collection  lens 
were  used  to  vary  the  maximum  magnon  wave  number  £max 
for  the  forward  scattered  light.  Second,  wave  vector  directions 
were  selected  though  the  use  of  a  I  mm  wide  slit-type  aperture 
placed  behind  the  collection  lens.  The  slit  orientation  relative 
to  the  applied  field  defines  the  usual  polar  spin  wave 
propagation  angle  6k . 

Figure  I  illustrates  the  scattering  geometry  and  the  wave- 
vector  selection  process  described  above.  The  Figure  shows  the 
YIG  film  propagation  structure  and  the  collection  optics.  The 
arrows  show  the  incident  light  with  wave  vector  q,  and  the 
scattered  light  at  q5  that  results  from  the  interaction  with  a 
magnon  with  wave  vector  kmax  .  This  upper  limit  on  |  kmax  | 
is  defined  by  the  radius  of  the  on-axis  diaphragm  shown 
schematically  on  the  far  right  side  of  the  diagram.  The  opening 
limits  the  scattered  light  into  the  spectrometer  to  those  photons 
that  derive  from  magnons  with  wave  numbers  ( k  )  less  than 
some  Amax  .  If  one  uses  the  1  mm  wide  slit  instead  of  the 
diaphragm,  as  indicated  by  the  dashed  lines,  then  one  can 
further  limit  the  selected  magnon  wave  vector  to  a  specific 
propagation  angle  6k  .  This  is  strictly  true  for  k  -values  above 
1000  rad/cm  or  so,  well  within  the  bounds  of  the  data 
considered  below.  A  full  slit  that  extends  over  the  entire 
aperture  of  the  collection  lens,  as  indicated  by  the  dashed 
rectangle,  was  used  for  the  initial  scans  to  determine  the  6k  - 
values  for  strong  scattering.  A  special  slit  that  extends  only  in 
one  direction  from  the  optic  axis  was  used  for  the  quantitative 
wave-vector  selective  measurements  to  be  considered  in  Sec. 
IV. 

For  the  wave-number  measurements  for  the  confluence 
magnons  in  section  IV,  an  additional  wave-vector  selection 
process  was  used.  In  this  set  up,  instead  of  a  diaphragm, 
different  size  on-axis  circular  light  blocks  were  used.  The 
diameter  of  the  block  plays  the  same  role  as  the  diameter  of  the 
diaphragm  in  Fig.  1,  except  than  one  now  selects  out  wave 
numbers  above  some  Amin  value  instead  of  k-  values  below 
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FIG.  1.  Schematic  diagram  of  the  wave-vector  selection  scattering 
geometry.  The  laser  light  is  incident  on  the  yttrium  iron  garnet  (YIG) 
Film  spin  wave  propagation  structure  at  wave  vector  q  ,  ,  with  a 
scattered  light  wave  vector  q  y  and  a  magnon  wave  vector  k  max 
oriented  at  an  angle  0max  relative  to  the  direction  of  the  applied  Field 
H  The  collection  lens  and  the  diaphragm/slit  assembly  define  the 
selected  magnon  wave  vector. 
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some  k  max  . 

Through  the  above  scheme,  one  can  perform  a  full  range  of 
wave  vector  selective  BLS  measurements.  Separate  data  runs 
for  different  diameter  (d )  diaphragms  or  circular  light  blocks, 
for  example,  can  allow  one  to  separate  out  the  scattering 
associated  with  specific  wave  numbers.  Diaphragms  with 

values  from  1  mm  to  30  mm  were  used,  typically  in 
increments  of  1  or  2  mm.  The  corresponding  £max  values  are 
1.2 xlO3  rad/ cm  and  3.6 xlO4  rad/cm,  respectively.  Light 
blocks  with  values  from  2  mm  to  10  mm  were  used.  These 
diameters  give  a  range  of  km]n  from  2.4x10^  to  1.2 xlO4 
rad/cm.  Measurements  were  always  made  with  the  same  light 
intensity  and  measuring  time.  Once  the  wave  number  has  been 
determined,  data  obtained  with  a  slit  aperture  can  be  used  to 
separate  out  the  product  magnons  as  a  function  of  propagation 
angle.  The  wave-vector  selectivity  capability  is  crucial  to  the 
identification  of  the  product  spin  waves  for  the  nonlinear  three 
magnon  splitting  and  confluence  processes  that  are  the  focus  of 
this  paper.  The  data  presented  in  Sec.  IV  will  illustrate  the  very 
different  wave-vector  make-up  for  the  parametrically  excited 
DESW  /p/2  magnons  from  the  splitting  process  and  the 
product  spin  waves  that  come  from  the  confluence  of  these 
magnons. 

Ill.  PARAMETRIC  DESW  AND  CONFLUENCE  MAGNONS 

How  do  the  BLS  data  reveal  the  presence  of  the  relevant 
magnons?  Figure  2  shows  representative  BLS  spectra  of 
scattering  intensity  vs.  frequency  for  different  input  cw  power 
levels,  as  indicated.  The  frequency  axis  shows  the  scattered 
photon  frequency  relative  to  the  laser  frequency.  The 
frequency  axis  goes  from  positive,  or  anti-Stokes  (AS),  on  the 
left,  to  negative,  or  Stokes  (S),  on  the  right.  The  strong 


Frequency  shift  [GHz] 

FIG  2.  Representative  Brillouin  light  scattering  spectra  of  scattering 
intensity  vs.  frequency  for  different  input  microwave  cw  power  levels 
(/jn  ),  as  indicated.  The  pump  frequency  (fp)  was  2.5  GHz.  The 
static  external  field  was  set  to  352  Oe.  The  nominal  spectrometer  free 
spectral  range  was  8  GHz..  There  were  no  wave-vector  selective 
diaphragms  in  place  and  the  measuring  time  was  the  same  for  all 
scans.  The  vertical  dashed  lines  and  labels  identify  different  peaks  in 
the  spectra  at  fp  ,  at  fc  ,  slightly  below  f p  ,  and  at  fpt 2. 


Rayleigh  peak  at  zero  frequency  comes  from  the  laser  line.  The 
incident  beam  was  positioned  in  the  middle  of  the  Y1G  strip 
about  2  mm  down  line  from  the  input  transducer.  The 
microwave  pumping  frequency  was  2.5  GHz.  The  static 
magnetic  field  was  set  at  352  Oe.  All  measurements  were 
made  with  the  same  incident  light  intensity,  the  wave  vector 
selective  diaphragm  removed,  and  the  same  TFP  scanning 
parameters. 

Note  also  that  the  scans  in  Fig.  2  show  magnon  signals 
only  on  the  positive  frequency  AS  side  of  the  spectra.  This  AS 
region  was  scanned  at  a  slow  scan  rate  in  order  to  bring  out  the 
magnon  structure.  Note  that  this  AS  region  corresponds  to 
scattering  in  which  one  of  the  microwave  or  parametrically 
pumped  magnon  is  destroyed  and  a  photon  is  created. 
Generally  speaking,  the  BLS  peaks  for  the  pumped  magnons  at 
fp  and  the  /p/2  split  magnons  are  found  to  be  stronger  on 
the  AS  side  of  the  spectra.  Intuitively,  this  is  consistent  with 
the  nature  of  the  scattering  from  magnons  driven  to  large 
occupation  numbers  by  either  the  microwave  pumping  process 
or  the  subsequent  interactions  that  produce  the  parametric 
magnons.  For  the  confluence  magnons,  however,  it  turns  out 
that  the  Stokes  as  well  the  Anti-Stokes  sides  of  the  spectra  are 
needed  to  obtain  a  full  picture  of  the  excited  modes.  Such  AS- 
S  confluence  magnon  effects  will  be  considered  in  Sec.  IV. 

The  spectra  examples  in  Fig.  2.  show  three  things.  (1) 
There  is  a  peak  in  the  spectra  at  the  pumping  frequency  that  is 
strong  at  low  power,  weakens  at  intermediate  power  levels,  and 
then  strengthens  again  at  the  highest  power  shown.  (2)  As  the 
power  is  increased,  one  sees  an  additional  peak  at 
/p/2«  1.25  GHz  appear  and  grow  in  strength.  (3)  One  can 
also  see  the  appearance  and  growth  of  a  third  peak  at 
fc « 2.25  GHz,  positioned  slightly  below  the  fp  peak  in 
frequency,  as  the  power  level  is  increased. 

In  reference  to  item  (1),  the  initial  drop  and  subsequent 
increase  in  intensity  for  the  peak  in  the  spectra  as  the  power 
level  is  increased  is  a  likely  indication  of  the  role  of  the  three 
magnon  processes  that  are  in  play  here.  Microwave  data  for 
manganese  ferrite  (MgF)  in  the  initial  Suhl  paper  on 
ferromagnetic  resonance  saturation  shows  a  gradual  decrease  in 
the  microwave  loss  as  one  moves  above  threshold  The  data  in 
Ref.  3,  however,  do  not  extend  to  the  extreme  25-  40  dB  above 
threshold  levels  needed  to  match  the  results  in  Fig.  2. 
Moreover,  the  nature  of  the  MgF  samples  is  not  specified.  One 
possible  scenario  is  that,  as  the  power  is  initially  increased  well 
above  threshold,  more  and  more  of  the  energy  of  the  pumped 
magnons  is  passed  to  the  parametric  magnons  and  the  peak  at 
/  drops  in  strength.  As  the  power  is  further  increased,  these 
parametric  processes  level  off  and  the  f  peak  grows  again. 

For  item  (2),  the  growth  of  the  fpj 2  magnon  splitting 
peak  with  increasing  power  provides  direct  evidence  for  the 
scenario  above.  This  is  actually  a  threshold  effect  similar  to 
that  found  for  Suhl  processes.  There  is  a  critical  power  level, 
typically  at  the  one  microwatt  level,  for  the  onset  of  the 
splitting  process.  Even  so,  the  splitting  signal  is  not  discernable 
from  the  BLS  data  until  the  input  power  is  much  higher,  as  in 
the  (b)  graph  of  Fig.  2.  As  elaborated  below,  these  product 
magnons  for  the  splitting  process  have  a  range  of  frequencies 
distributed  around  fp  jl .  The  growth  in  the  peak  at  fc  [item 
(3)]  slightly  below  fp  is  ascribed  to  the  confluence  of  some 
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subset  of  these  split  magnons.  The  concluding  discussion  for 
this  section,  given  below,  provides  qualitative  connections 
between  these  data  and  the  proposed  magnon  splitting  and 
confluence  processes.  Quantitative  connections  will  follow  in 
Sec.  IV,  based  on  detailed  wave-vector  selective  measurements. 

These  connections  are  based  on  the  schematic  MSBVW 
spin  wave  dispersion  diagrams  and  specific  spin-wave  mode 
points  shown  in  Fig.  3.  A  detailed  description  of  the 
nomenclature  is  given  in  the  expanded  Figure  caption.  The 
curves  were  computed  from  the  standard  DESW  theory  in  Ref. 
16,  standard  YIG  parameters,  and  the  same  film  thickness  and 
static  field  cited  for  the  measurements.  The  identified  mode 
points  are  intended  to  make  qualitative  connections  with  the 
BLS  results. 

Note  that  the  actual  parametric  modes  need  not  have  wave 
vectors  that  are  strictly  co-linear  with  the  field.  Keep  in  mind 
that  a  given  MSBVW  pumped  magnon  at  /  and  wave  vector 
k  will  split  into  two  parametric  magnons  with  frequencies 
fs[  and  fs 2  and  wave  vectors  ksl  and  kj2  i  respectively. 
Energy  and  momentum  conservation  impose  the  conditions 
fp  =  /,  1  +fs2  and  kp  =  +  k,2  •  since  \  kp  |  is  very  small 

compared  to  |  k5l  |  and  |  ksl  | ,  it  is  clear  that  k5l  and  k52  must 
be  nearly  equal  and  oppositely  directed.  The  dispersion  curve 
constraints  also  force  both  ksl  and  kj2  to  be  nearly  co-linear 
with  the  field.  Similar  remarks  apply  to  the  reverse  confluence 
process  that  take  various  pairs  of  split  magnons  into  modes  at 
frequencies  close  to  fc .  The  up  shift  in  frequency  to  fk  «  fc 
in  this  case  leads  to  allowed  modes  at  relatively  large 
propagation  angles.  As  the  data  of  the  next  section  will  show, 
the  parametric  spin  waves  turn  out  to  have  in-plane  propagation 
angles  away  from  the  field  direction  by  ±10°  or  so  for  the  split 
fpj 2  magnons  and  by  ±45°  or  so  for  the  confluence  fc 
magnons. 

Overall,  Fig.  3  shows  a  consistent  map  of  the  connections 
between  the  initially  excited  cw  MSBVW  spin  wave  signal  in 
the  YIG  strip  and  the  parametric  magnons  observed  from  the 
BLS  measurements.  If  one  considers  a  low  wave-number 
MSBVW  mode  at  fp  (upper  open  diamond),  magnon  splitting 
processes  that  conserve  frequency  and  momentum  give  product 
magnons  as  indicated  by  the  square  and  circle  points  clustered 
around  the  fp! 2  line  cut  in  the  figure.  At  the  same  time,  the 
lower  frequency  subset  of  these  split  magnons  can  be  combined 
to  give  a  confluence  peak  at  fc  .  The  dashed  box  in  the  lower 
central  part  of  Fig.  3  identifies  two  such  modes. 

The  layout  of  modes  in  Fig.  3  leads  to  one  more  important 
connection  with  the  data.  Note  that  both  the  splitting  peak  in 
Fig.  2  at  fpj 2  and  the  confluence  peak  at  fc  are  rather  broad. 
Based  on  the  scenario  given  above,  one  can  explain  the  f  jl 
line  broadening  in  terms  of  the  relatively  wide  band  of  split 
magnons  shown  by  the  circle  and  square  points  in  Fig.  3.  It 
also  stands  to  reason  that  if  the  observed  split  magnon  peak  is 
broad,  any  peak  due  to  the  confluence  of  these  magnons  will 
also  be  broad. 

Taken  at  face  value,  the  connections  implied  in  Fig.  3  are 
nothing  more  than  conjecture.  The  mode  connections  are  based 
entirely  on  the  dispersion  curves  shown  and  the  frequencies 
match  those  found  in  the  experiment.  What  remains  is  to  show, 
from  actual  wave-vector  selective  BLS  data,  that  the  split 
modes  clustered  around  fpjl  in  frequency  and  the  product  low 


FIG.  3.  Schematic  dispersion  curves  of  spin  wave  frequency  fk  vs. 
wave  number  k  with  mode  connections.  The  curves  are  for  the  first 
six  dipole-exchange  spin  wave  branches  for  in-plane  propagation 
parallel  or  anti-parallel  to  the  field  H  and  an  H  -value  of  352  Oe,  the 
same  as  for  the  data  in  Fig.  2.  The  horizontal  dashed  lines  label  the 
frequencies  for  the  three  peaks  identified  in  Fig  2  for  the  directly 
pumped  magnetostatic  backward  volume  wave  (MSBVW)  spin  waves 
at  fp  =  2.5  GHz,  confluence  magnons  at  fc  ,  and  splitting  magnons 
at  fp! 2  .  The  open  diamond  at  fp  labels  the  MSBVW  operating 
point.  The  matching  pairs  of  solid  and  open  circles  and  squares  show, 
in  schematic  form,  parametrically  excited  split  magnons  at  ±k  and 
with  fk  values  equally  shifted  up  and  down  from  the  /  /2  line.  The 
open  diamonds  for  a  range  of  low  |  k  |  values  on  the  fc  line  show  the 
nominal  position  in  frequency  and  wave  number  for  the  confluence 
magnons.  The  dashed  square  close  to  1.12  GHz  identifies  split 
magnons  that  give  a  confluence  magnon  at  fc  . 

wave-number  confluence  modes  at  fc  match  the  data. 

IV.  WAVE  VECTOR  SELECTIVE  ANALYSIS 

The  techniques  and  procedures  for  the  wave-vector 
selective  BLS  analysis  were  outlined  in  Sec.  II.  Recall  that 
there  are  two  stages  in  the  analysis  for  the  split  magnons,  one 
with  a  variable  diameter  diaphragm  for  wave-number  analysis 
and  one  with  a  rotating  slit  aperture  for  wave-vector  direction 
analysis.  Recall  the  key  feature  of  the  stage  one  measurements, 
namely,  the  clear  and  distinct  fpjl  peak  in  Fig.  2(d)  that 
appears  suddenly  when  the  diameter  of  the  wave-vector 
selective  aperture  is  larger  than  about  9  mm.  This  critical 
d-  value  corresponds  to  &max  «  1  .lx  104  rad/cm  .  This  is 
precisely  the  minimum  allowed  /revalue  for  f  j 2  magnons 
from  DESW  theory.  This  match-up  comprises  prima  face 
evidence  for  the  connection  with  DESW  magnons  as  proposed 
above. 

Figure  4  shows  representative  wave-vector  selective  BLS 
data  for  the  /  jl  peak  in  Fig.  2  and  companion  iso-frequency 
magnon  wave-vector  contour  plots  for  the  low  order  DESW 
spin-wave  dispersion  branches  at  f k  -  /  /2  .  Graph  (a)  shows 
the  BLS  signal  as  a  function  of  the  cut-off  &max  wave  number 
value.  The  vertical  dashed  lines  denote  the  upper  and  lower 
limit  A:max  values  for  the  diaphragm  set-up  described  above. 
Figure  4(b)  shows  calculated  iso-frequency  contour  plots  for 
the  1.25 GHz  cut  of  the  dipole-exchange  spin  wave  dispersion 
surface.  While  the  horizontal  axes  for  both  graphs  correspond 
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FIG.  4.  (a)  BLS  scattering  signal  at  fpj2  as  a  function  of  the 
maximum  allowed  wave  number  £max  for  an  input  power  of 
5.56mW  and  fp-2.5QHz.  The  vertical  dashed  lines  denote  the 
limit  kmax  values  for  the  experimental  diaphragm  set-up.  (b) 
Theoretical  iso-frequency  contour  plots  of  the  transverse  in-plane  spin 
wave  wave-vector  component  /rj_  vs.  the  parallel  component  k //  for 
the  four  lowest  order  DESW  modes  at  fpj 2  .  The  plot  was  obtained 
from  the  DESW  theory  of  Ref  16,  with  parameters  that  match  the 
experiment.  The  open  circles  and  k  q  labels  mark  the  low  limit  for  the 
available  DESW  spin  waves  at  fp /2  .  The  angle  labeled  in  (b) 
shows  the  maximum  allowed  spin  wave  propagation  angle 

to  wave  number,  the  plotted  parameter  in  (b)  corresponds  to 
the  component  of  the  in-plane  DESW  wave  vector  k  that  is 
parallel  to  the  static  field  H  ,  labeled  as  k), .  For  comparison 
purposes,  the  horizontal  scales  for  (a)  and  (b)  are  the  same. 
The  vertical  scale  in  (b)  corresponds  to  the  in-plane  transverse 
component  of  k  ,  labeled  as  .  For  reference  purposes  for  the 
discussion  of  the  angle  data  to  follow,  the  0™™  labeled  in  (b) 
indicates  the  maximum  theoretical  in-plane  spin-wave 
propagation  angle  for  the  half-frequency  spin  waves. 

The  data  in  (a)  and  the  iso-frequency  contours  in  (b) 
confirm  the  match-up  noted  in  the  introduction  to  this  section. 
The  onset  of  detected  spin  waves  as  a  function  of  the 
experimental  /tmax  value  occurs  precisely  at  the  theoretical 
lower  limit  k  value  at  k0  «  1 .1  x  104  rad/cm  for  the  half¬ 
frequency  DESWr  spin  waves.  These  lower  limit  k-  values  are 
marked  in  (b)  by  open  circles.  Note  that  this  onset  point  in  (b) 
is  for  a  DESW  propagation  direction  that  is  parallel  to  the  static 
field  and  with  a  propagation  angle  0k  equal  to  zero  or  180° 
only.  The  full  iso-frequency  contours  in  (b),  however,  provide 
much  more  information  that  the  confirmation  of  this  onset  point 
for  split  magnons.  The  contours  allow  one  to  read  off  the 
allowed  propagation  angles  for  the  split  magnons  with 
k>k0  =  1.1  x  104  rad/cm.  From  these  plots,  one  expects  that  the 
angle  selective  BLS  data  will  reveal  spin-wave  propagation 
angles  for  the  half-frequency  spin  waves  in  the  range 
-0™  <  ek  <0™  and  (1 80°  -  0™™  )<  @k  <(180°  +  <9™)- 

As  noted  previously,  0k  resolved  BLS  data  were  obtained 
with  a  rotating  slit-type  aperture  scheme.  In  order  to  resolve 
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FIG.  5.  Polar  representation  of  the  BLS  signal  for  the  fpj  2  DESW 
magnon  peak  as  a  function  of  the  slit  aperture  rotation  angle  9k  .  The 
arrow  shows  the  direction  of  the  static  field  H  .  The  dashed  lines 
denote  the  double  sector  of  spin  wave  propagation  angles  for 
and  (1 80°  -  0™ )  <  0*  <  (1 80°  +  0™ )  that  match 
the  range  of  allowed  9k  values  from  Fig  4b. 

the  full  0-360°  range  of  possible  0  k  —  values,  the 
measurements  were  done  with  a  half  slit  arrangement  that 
extends  outward  from  the  optical  axis  in  one  direction  only, 
rather  than  the  full  slit  shown  in  Fig.  1.  Figure  5  shows  a 
representative  polar  plot  of  BLS  intensity  data  as  a  function  of 
0k  for  the  /p/2  split  magnon  peak  in  Fig.  2(d)  at 
Pm  =  5.56  mW.  As  in  Fig.  4,  the  field  is  parallel  to  the  0k  =  0 
direction.  The  data  were  obtained  for  angular  steps  of  2  and 
with  the  same  scanning  times  for  all  measurements.  The  dashed 
lines  at  angles  of  0j™  to  the  horizontal  mark  the  expected 
range  of  pumped  magnons  from  Fig.  4.  The  angular 
distribution  of  the  DESW  /  / 2  magnons  is  clear  from  the 
figure  and  no  extensive  discussion  is  needed  here.  The  data 
show,  rather  unambiguously,  that  the  BLS  signal  is  constrained 
within  a  double  sector  of  angles  within  ±0™  of  the  field  axis. 
The  span  of  observed  propagation  angles  matches  nicely  to  the 
contour  plot  in  Fig  4(b). 

The  above  data  and  theoretical  connections  provide  a 
complete  map  of  the  DESW  fpj 2  magnons  that  result  from 
the  three  wave  splitting  process  for  pumped  MSBVW  magnons. 
Both  the  directly  measured  wave  numbers  and  the  propagation 
directions  of  the  split  magnons  match  nicely  to  theory.  This 
represents  the  first  direct  experimental  map  of  three  magnon 
splitting  processes  in  any  magnetic  system  pumped  in  the 
backward  wave  configuration. 

Turn  now  to  the  confluence  magnons.  Recall  that  in  this 
case,  blocking  circles  as  well  as  diaphragms  and  slit  apertures 
were  used  for  the  wave-vector  selective  measurements.  The 
changes  in  intensity  with  diaphragm  diameter  d  for  the  fc 
peak  are  quite  different  from  those  for  the  split  magnons.  In 
this  case,  the  signal  increases  steadily  as  d  is  increased  from  I 
to  6  mm  and  then  appears  to  increase  at  a  slightly  lower  rate  up 
to  the  limit  d  -value  of  10  mm.  This  implies  that  the 
predominant  confluence  magnon  k-  values  are  generally  below 
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about  IO'racL'cm  or  so.  For  the  light  block  in  place  of  the 
diaphragm,  the  fc  peak  signal  intensity  starts  high  and 
gradually  diminishes  to  near  zero  as  d  is  increased  from  the 
initial  2  mm  value  to  the  10  mm  upper  limit  value.  This  serves 
to  confirm  a  confluence  magnon  distribution  with 
/c  <  1 04  rad/cm  and  predominant  k  -values  somewhat  below 
this  limit.  At  the  same  time,  the  slit  aperture  data  indicate 
propagation  angles  that  are  typically  30°-60°  and  I20°-I50° 
away  from  the  field  axis.  Recall  that  the  split  magnons  are  all 
very  close  to  the  field  axis.  It  is  clear  from  these  results  that  the 
wave  vector  make-up  of  the  confluence  magnons  is  very 
different  from  that  for  the  split  magnons. 

Figure  6  shows  representative  wave-vector  selective  BLS 
data  for  the  fc  peak  in  Fig.  2  and  companion  iso-frequency 
magnon  wave  vector  contour  plots  for  the  low  order  DESW 
spin-wave  dispersion  branches  at  fk=fc-  Graph  (a)  shows  the 
fc  peak  intensity  as  a  function  of  the  diaphragm  and  light 
block  cut-off  wave  numbers,  £max  and  kmm  ,  as  indicated. 
The  shaded  curve  at  the  bottom  of  the  graph  provides  a 
qualitative  indication  of  the  fc  magnon  k  -distribution  implied 
by  the  two  sets  of  data.  Graph  (b)  shows  calculated  iso¬ 
frequency  contour  plots  for  the  2.25  GHz  fc  cut  of  the  dipole- 
exchange  spin-wave  dispersion  surface.  The  angle  0J^ 
indicates  the  maximum  theoretical  in-plane  spin  wave 
propagation  angle  for  the  confluence  magnons.  Graph  (c) 
shows  representative  BLS  data  in  a  polar  plot  representation  for 
the  fc  confluence  magnon  peak  in  Fig.  2.  In  order  to  obtain 
the  somewhat  balanced  4-lobe  plots  shown  in  (c),  it  was 
necessary  to  use  both  anti-Stokes  positive  frequency  side  BLS 
data  for  90°  <  Gk  <  270°  and  Stokes  negative  frequency  side 
data  for  ~90°(270°)  <  9k  <  90°  . 

The  results  in  Fig.  6  provide  a  useful  map  of  the  magnon 
make-up  for  the  BLS  confluence  peak.  These  results  also 
demonstrate  the  marked  difference  between  these  confluence 
magnons  at  fc  and  the  split  magnons  at  fp! 2  discussed  at  the 
start  of  this  section.  Graph  (a),  for  example,  shows  that  the  fc 
magnons  are  limited  to  k  values  below  IO4  rad/cm  ,  while  the 


data  in  Fig.  4  show  that  the  split  magnons  are  all  at  wave 
numbers  well  in  excess  of  1 04  rad/cm  .  Graph  (b)  shows, 
moreover,  that  the  expected  angular  distribution  of  the  fc 
magnons  emphasizes  propagation  directions  at  about  45°  or  so, 
relative  to  the  field  axis,  while  the  predominant  0k  values  for 
the  split  magnons  are  much  closer  to  the  field  axis.  The  angular 
data  in  (c)  confirm  that  the  magnons  associated  with  the  fc 
peak  do,  in  fact,  give  polar  plot  lobes  that  match  the  theory. 

V.  PUMPING  FREQUENCY  DEPENDENCES 

The  results  in  Secs.  Ill  and  IV  serve  to  quantify  the 
presence  of  (1)  DESW  fpj 2  magnons  that  result  from  the 
parametric  three  wave  splitting  process  for  pumped  MSBVW 
magnons,  and  (2)  confluence  magnons  at  f  that  result  from 
the  recombination  of  the  DESW  magnons.  The  wave-number 
make-up  and  propagation  characteristics  of  these  magnons  were 
mapped  experimentally  and  connected  to  theory.  All  of  these 
results,  however,  are  based  on  data  for  one  fixed  value  of  the 
in-plane  field  and  one  specific  drive  frequency  selected  for 
strong  coupling  to  the  pumped  MSBVW  mode.  This  section  is 
concerned  with  one  further  test  of  the  mode  connections 
established  above,  based  on  the  effect  of  changes  in  the 
pumping  frequency  f  on  the  BLS  response,  with  the  field 
held  at  the  same  value  as  before.  The  strategy  here  is  similar  to 
that  of  Synogach  et  al.  in  Ref.  1.  In  that  work,  done  for  the 
case  of  pumped  magnetostatic  surface  modes,  an  increase  in 
both  the  field  and  the  pumping  frequency  resulted  in  an  upward 
shift  in  both  the  spin  wave  band  and  the  operating  point 
frequency  to  the  point  where  three  magnon  splitting  processes 
are  no  longer  allowed.  This  elimination  of  three  magnon 
splitting  processes  resulted  in  a  major  change  in  the  microwave 
response  of  the  surface  wave  delay  line  structure. 

The  approach  here  is  much  simpler.  The  field  was  held  at 
the  same  value  as  before,  but  now,  the  pumping  frequency  was 
gradually  decreased  to  the  point  where  the  downshifted  fpl  2 
cut  across  the  spin  wave  dispersion  curve  diagram  falls  below 
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FIG.  6.  (a)  Data  on  the  BLS  scattering  signal  level  at  the  fc  peak  vs.  the  two  limit  wave  numbers  A:min  and  kmax  ,  as  obtained  for  different  diameter 
light  blocks  or  diaphragms,  respectively.  The  dashed  light  gray  area  gives  a  qualitative  estimate  of  the  confluence  magnon  k  -  distribution  implied  by 
these  data  (b)  Theoretical  iso-frequency  contour  plots  of  the  transverse  in-plane  spin  wave  wave-vector  component  £j_  vs.  the  parallel  component 
k  jf ,  as  referenced  to  the  field  direction,  for  the  five  lowest  order  DESW  modes  at  fc  .  The  plot  was  obtained  from  the  DESW  theory  of  Ref  16,  with 
parameters  that  match  the  experiment.  The  angle  9™^*  corresponds  to  the  maximum  allowed  spin-wave  propagation  angle  for  confluence  magnons  at 
the  fc  frequency  cut.  (c)  Polar  representation  of  the  BLS  signal  level  at  the  fc  peak  as  a  function  of  the  slit  aperture  rotation  angle  for  the  input 
power  of  5.56  mW.  The  arrow  shows  the  direction  of  the  static  magnetic  bias  field  H  . 
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the  bottom  of  the  band  of  DESW  modes.  One  can  see  this 
basic  effect  from  Fig.  3.  The  bottom  of  the  band  sits  at  about 
1.0125  GHz.  This  means  that  a  pumping  frequency  below 
about  twice  this  value,  or  2.025  GHz,  will  be  low  enough  to 
eliminate  the  possibility  of  any  three  magnon  splitting  response 
at  all.  One  can  also  see  that  as  the  pumping  frequency  is 
gradually  lowered,  the  fp/2  cut  of  modes  across  the  bottom  of 
the  band  will  become  narrower  and  narrower  and  the  mean 
k—  value  for  the  split  magnons  will  shift  out  to  about 
4xl04  rad/cm  just  before  the  splitting  response  vanishes.  This 
will  also  cause  the  fc  point  for  the  confluence  magnons  to 
move  closer  and  closer  to  fp  . 

Panels  (a)  -  (d)  in  Fig.  7  show  the  effect  of  a  lowering  of 
the  pump  frequency  on  the  spectra  in  Fig.  2.  The  data  go  from 
the  range  where  three  magnon  splitting  is  readily  allowed,  as  in 
(a)  for  fp  =  2.4  GHz,  to  a  range  where  such  processes  are 
barely  allowed,  as  in  (d)  for  fp-2A  GHz.  All  data  were 
obtained  for  the  same  configuration  and  parameters  as  for  the 
data  in  Fig.  2(d),  except  for  the  shift  in  pumping  frequency. 
The  three  vertical  dashed  lines  in  (a)  and  (b)  mark  the  pumping 
frequency  operating  point  ( fp )  ,  the  confluence  (/c)  and 
splitting  (f  1 2)  peaks  from  the  data,  as  indicated.  In  (c),  the 
two  vertical  dashed  lines  mark  the  f  and  /  / 2  points  only. 
In  this  case,  the  fc  point  is  so  close  to  the  operating  point  that 
the  two  frequencies  are  indistinguishable.  The  single  vertical 
dashed  line  in  (d)  marks  the  f  point  only.  The  arbitrary  unit 
[a.u.]  vertical  axis  BLS  intensity  scales  for  the  four  graphs 
show  self  consistent  relative  values.  The  matching  scales  for 
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FIG.  7.  Representative  BLS  spectra  of  scattering  intensity  vs. 
frequency  for  different  pumping  frequencies,  (a)  f  =  2.4  GHz ,  (b) 
fp  =  2.3  GHz  ,  (c)  fp  =  2.2  GHz  and  (d)  fp  =  2.1  GHz  .  The  static 
external  field  was  set  to  352  Oe.  The  nominal  spectrometer  free 
spectral  range  was  8  GHz.  There  were  no  wave-vector  selective 
diaphragms  or  light  blocks  in  place  and  the  measuring  time  was  the 
same  for  all  scans.  The  vertical  dashed  lines  and  labels  identify 
different  peaks  in  the  spectra  at  fp  ,  at  fc  ,  and  at  fD! 2. 


(a),  (b),  and  (c)  demonstrate  the  evolution  of  the  spectra  and  the 
individual  intensities  with  the  drop  in  frequency.  The  larger 
range  for  (d)  shows  the  rapid  rise  in  the  strength  of  the  fp  peak 
as  one  approaches  the  cut-off  point  for  splitting. 

The  four  panels  in  Fig.  7  demonstrate  the  evolution 
scenario  outlined  above.  Recall  that  as  fp  is  gradually  reduced, 
the  gap  between  the  fpt  2  and  the  bottom  of  the  band  becomes 
smaller  and  smaller.  This  means  that  the  fc  point  must  also 
move  closer  and  closer  to  fp  .  This  is  precisely  what  one  sees 
from  panels  (a)  -  (c)  in  Fig.  7.  For  (a)  the  fp  to  fc  spacing  has 
dropped  from  the  nominal  0.25  GHz  value  in  Fig.  2  to  0.1  GHz. 
For  (b)  this  spacing  has  dropped  further  to  about  0.05  GHz.  In 
(c),  the  fc  peak  position  is  not  even  discemable  from  fp  on 
the  scale  of  the  graph.  Finally,  for  the  case  shown  in  Fig.  7(d), 
with  the  bottom  of  the  band  at  about  1.02  GHz  and  /_/ 2  at 
1.05  GHz,  the  half  frequency  point  is  so  close  to  the  bottom  of 
the  band  that  the  splitting  and  confluence  effects  of  interest 
here  no  longer  occur.  In  panel  (d),  one  sees  that  the  fpl2 
splitting  peak  has  disappeared  entirely  and  the  peak  at  fp  has 
increased  in  intensity  by  a  factor  of  3  or  so. 

The  above  points  mainly  relate  to  the  evolution  in  the 
shape  of  the  spectra  as  the  pumping  frequency  is  decreased. 
Note  also  that  the  data  show  a  clear  degradation  in  the  shape  of 
the  fpl2  splitting  peak  as  one  goes  from  panel  (a)  to  panel  (c). 
This  degradation  can  be  attributed,  quite  plausibly,  to  two 
things,  (1)  the  move  of  the  fpl 2  frequency  cut  closer  and 
closer  to  the  bottom  of  the  band,  and  (2)  the  wave  number 
limits  in  the  collection  optics.  From  Fig.  3,  one  can  see  that  as 
one  moves  from  the  situation  shown  to  a  pump  frequency  at  the 
bottom  of  the  band,  the  k  -values  for  the  split  frequency 
modes  that  contribute  to  the  confluence  will  move  from  the 
indicated  value  of  about  2xl04  rad/cm  out  to  about 
4xl04  rad/cm  .  At  the  same  time,  the  maximum  A: -value  that 
can  be  accessed  by  the  collection  lens  is  about  3.6 xl(T  rad/cm  . 
This  means  that  the  somewhat  washed  out  fpl2  peaks  in  (b) 
and  (c)  are  for  wave  numbers  that  approach  or  exceed  the  limits 
of  the  set-up. 


V.  SUMMARY  AND  CONCLUSION 

Three  magnon  splitting  and  confluence  processes  in  thin 
YIG  films  for  the  MSBVW  propagation  configuration  has  been 
experimentally  observed  by  Brillouin  light  scattering 
measurements  and  confirmed  by  wave-vector  selective  analysis. 
A  complete  characterization  for  both  of  these  nonlinear 
processes  has  been  made.  The  change  in  the  response  with 
pumping  frequency  as  well  as  the  fully  analyzed  wave-vector 
make  up  of  the  product  magnons  provide  proof  positive  of 
these  processes.  Interestingly,  and  perhaps  fortuitously,  the 
data  reveal  a  confluence  response  for  only  selected  half¬ 
frequency  split  magnons  that  are  slightly  downshifted  from  the 
strict  half-frequency  point.  This  state  of  affairs  pushes  the 
confluence  peak  to  a  point  slightly  below  the  peak  at  the  signal 
frequency  and  renders  it  observable. 
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